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Abstract 

Bosonization of degenerate fermions yields insight both into Landau Fermi liquids, 
and into non-Fermi liquids. We begin our review with a pedagogical introduction 
to bosonization, emphasizing its applicability in spatial dimensions greater than one. 
After a brief historical overview, we present the essentials of the method. Well known 
results of Landau theory are recovered, demonstrating that this new tool of many-body 
theory is robust. Limits of multidimensional bosonization are tested by considering 
several examples of non-Fermi liquids, in particular the composite fermion theory of 
the half-filled Landau level. Nested Fermi surfaces present a different challenge, and 
these may be relevant in the cuprate superconductors. We conclude by discussing the 
future of multidimensional bosonization. [Published in Advances in Physics 49, 141 - 
228 (2000).] 
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1 Introduction 

Bosonization, a concept which sounds quite mysterious, is in fact both easy to understand 
and ubiquitous. In this review we outline an approach via bosonization to the problem of 
many interacting fermions in spatial dimensions greater than one. As the virtue of multi- 
dimensional bosonization lies in the reduction of complicated four-fermion interactions to a 
Gaussian problem, we stress the simplicity of the approach. 

Multidimensional bosonization is more delicate than bosonization in one dimension. In 
particular the technique is not useful for certain strongly correlated systems, specifically 
those in which the Fermi surface is obliterated by singular or unscreened interactions, as 
it relies heavily on the existence of an abrupt change in the quasiparticle occupancy near 
the Fermi wavevector. It is also difficult to implement systematic perturbative expansions 
using bosonization. In these instances, it is more useful to work with the original fermionic 
variables. Nevertheless bosonization provides a straightforward approach to many interesting 
problems and we begin this review with a brief introduction to the method in its several 
forms. 

Historically bosonization has been associated closely with the quantum physics of parti- 
cles moving in one spatial dimension or with statistical mechanics in two dimensions. (For 
up-to-date reviews of one-dimensional bosonization, which contain material not covered here, 
see the articles by Schulz[|l|, Voit0, and von Delft and Schoeller||.) But any fermion system 
which exhibits collective excitations has some sort of bosonic description. Spin waves in a 
ferro magnet, plasmons in a metal, and superconducting order are well known examples. In 
these instances, some if not all of the low energy degrees of freedom have a bosonic descrip- 
tion. We focus here on cases when all of the low-energy fermion degrees of freedom may 
be replaced by bosons. Mathematically the replacement can be carried out, at least for- 
mally, by means of a Hubbard- Stratonovich transformation^: auxiliary bosons are coupled 
to the fermion degrees of freedom, and then the fermions are integrated out yielding a purely 
bosonic effective theory. In practice, however, the transformation is useful only if there are 
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additional constraints. The following two examples illustrate the power of bosonization when 
the fermions are subject to constraints in the form of infinite symmetries. 

Consider first the problem of representing a spin-1/2 degree of freedom in terms of the 
underlying fermions. The angular momentum algebra is reproduced faithfully if we introduce 
fermion creation and annihilation operators c' a and c a which obey the canonical anticom- 
mutation relations {c^ a , Cg} = 8p. Then we may write: 

S = \c^ a c p (1.1) 

subject to the constraint 

n = c ]a c a = 1 . (1.2) 

Here an implicit summation over repeated raised and lower Greek spin indices is assumed. 
More generally, to describe a quantum magnet for instance, we attach a site index x to the 
operators: 

5x = ~4« ai c x/3 . (1.3) 
If, for example, we consider now the nearest-neighbor Heisenberg magnet 

H = J £ 4-4 (1-4) 

<x,y> 

it is clear that its representation in terms of the underlying fermions is not as economical as 



its representation in terms of the spins themselves, since the spin operators of equation |L3 



and hence the Hamiltonian equation O, are invariant under U(l) gauge transformations at 
each lattice point |5], ||: 

r _> J 8 * c 

4° - e-^cf (1.5) 

because the local phase rotation # x cancels out. We denote this infinite gauge symmetry 
£7(1)^°. The physical meaning of this local symmetry is simply that the charge degrees of 
freedom are frozen out and play no role in the insulating magnet. The underlying fermions 
have both charge and spin degrees of freedom, but as the number of fermions on each site is 
fixed to be one, only the spin degree of freedom is active. 

Spin operators S£ and S~ create particle-hole excitations, localized on a single site x. 
That is, S x creates a spin-j fermion and destroys a spin-j fermion at site x. These local 
particle-hole excitations, taken together, have a bosonic character, as two fermions bound 
together obey boson statistics. As the charge degree of freedom is frozen in place, the 
remaining spin degree of freedom can be represented instead in terms of bosons which obey 
the usual canonical commutation relations [b a , b^} = Now it is permissible to write: 

& = ^ a 32 I*. (1-6) 

These spins obey precisely the same algebra, and if we now constrain the number of bosons 
to be one per site, 

n x = fc x a & xa = 1 , (1.7) 
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the operators are spin-|. Note that such a replacement of fermions with bosons - our first 
concrete example of bosonization - would not have been permissible for a model which 
contains both charge and spin degrees of freedom, such as the Hubbard model. For instance, 
in the non- interacting U — > limit, the Hubbard model reduces to the tight-binding model 

H = -t Y, ci a c ya + H.c. (1.8) 

<x,y> 

which is clearly not invariant under the U(l)™ transformations, equation |1.5| . Distinct U(l) 
rotation angles at neighboring sites, 8 X and 9 y do not cancel in general. Physically this means 
that the full Fermi character of the electrons is manifested, as they propagate from lattice 
site to lattice site. Of course the ground state in the U — > limit is a degenerate Fermi 
liquid. For U > at half-filling, and for a bipartite lattice, on the other hand, the charge 
degrees of freedom are frozen out. Then the low energy effective theory of the resulting 
antiferromagnetic insulator does exhibit the ^(1)^° symmetry. 

At this point it might be asked what practical benefit can be gained by the replacement of 
the fermions with bosons. Our model provides a nice illustration: with bosons, semiclassical 
approximations may be contemplated. In the case of the antiferromagnet with J > 0, for 
instance, an expansion in powers of 1/S about the limit of large on-site spin S is now possible 
jTj, as the on-site spin can be increased simply by populating more bosons at each site: 
S = \rib- (In the fermion representation this can only be done by attaching an additional 
orbital index to the fermions, and then enforcing Hund's rule to align the fermions into a 
maximal spin state.) In particular, the Neel ordered state of the unfrustrated antiferromagnet 
on a square lattice can be accessed in this way even at the mean-field level. In contrast mean- 
field theories with fermions are disordered 0. 

There is another system which possesses infinite U(l) symmetry and which can be 
bosonized easily: the free, degenerate, Fermi gas or liquid [§]. We begin with the usual 
Hamiltonian: 

^ = E e kcj c a c ka . (1.9) 

k 

The Hamiltonian is invariant under an infinity of U(l) rotations, but now at each point in 
momentum space, rather than position space: 



Cka * 6 Ckct ' 



1.10) 



This U(l)^ symmetry, like the U(l)^ symmetry explored above, has a clear physical in- 
terpretation: in the absence of interactions and disorder, momentum is a good quantum 
number for the single-particle states, as these are infinitely long-lived. Thus we may make 
the identification: 

£/(l)^° <-> insulating solid 

U(l)™ <-> conducting liquid. (1-11) 

Again the particle-hole pairs can be bosonized, and this bosonization furnishes us with 
another example of how some properties can be seen more clearly in the bosonic basis. Recall 
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the classic calculation of the specific heat at constant volume V, Cy, of a degenerate Fermi 
liquid||. Two integrals, each a function of the temperature T and the chemical potential /z, 
must be computed: the total energy of the system U (T, //) and the total number of fermions 
N(T,fi). To extract CV, the chemical potential must be determined first, as it depends on 
temperature. Taking this variation into account, the correct answer at low temperatures is: 



C v =(%) =^*f,IW(0)xV. (1.12) 



Here N(0) is the density of states per unit volume at the Fermi surface. The same result 
may be recovered more readily by simply treating the particle-hole excitations as free bosons. 
More readily, because there is no chemical potential for the bosons as they are not conserved 
in number, but rather vanish in the zero-temperature limit. Now the total energy of a system 
of free bosons with single-particle energies e« is given by: 

u = y—^ — (i.i3) 

t-J e ti/k B T _ I V > 

and therefore 

C " = 4^?^5)' (L14) 

To apply this formula we need a sensible interpretation of e^. As we are interested in the 
low-temperature properties, we first linearize the spectrum of particle-hole excitations: 

(k + q) 2 -k 2 



2m 

k • q 



m 



;i.i5) 



where the second line holds in the limit |q| <C |k|. This restriction means that we are consid- 
ering only particle-hole excitations where particle and hole have nearly the same momentum. 
We show below that this constraint, which is unphysical, can in fact be discarded. At low 
temperatures the particle-hole pairs lie near the Fermi surface, and we may set |k| = kp, the 
Fermi momentum. Then the dispersion reads 

ek+q - e k = v F q\\ (1.16) 

where q» is the component of the momentum parallel to the Fermi wavevector k at the Fermi 
surface, q» = k • q. Of particular significance is the fact that the energy of a pair has this 
simple form when |q| <C hp, not depending on k, but rather only on the relative momentum 
q» . The energy of a particle-hole pair, like that of a single free particle, is determined uniquely 
by a single momentum. Substituting the dispersion into the specific heat formula, equation 
|1.14j , we obtain: 

No. of states on the Fermi Surface ^ (f^gy) 2 
Cv = JTZ^i x 



^bT 2 ^osinh 2 (30 



8nk F (L/2n) 3 [°° ^ v F q\ 



Ak B T> - Jo sinh 2 ^) • (L17) 



5 



Upon making a change of variables, and using the fact that 

x 2 dx IT 2 



1.18) 



h sinh 2 (x) 6 

the exact result, equation |1.12| , is recovered with the correct prefactorp|. 

Why is it permissible to restrict |q| <C kpl At any non-zero temperature there are a 
macroscopic number of particle-hole excitations, some with large q and some with small q. 
However, in the vicinity of any point on the Fermi surface there are statistically equal num- 
bers of particles and holes. Whether or not these particles and holes originate from small or 
large momentum scattering processes, they have the same statistical mechanical description, 
and they can always be re-assigned to small-q pairs which are described accurately by the 
above bosonic calculation. 

With this introduction to multidimensional bosonization, we turn now to a brief history of 
efforts to use bosonization to describe interacting fermion liquids. For more than forty years 
Landau Fermi liquid theory provided the framework to discuss the physics of strongly inter- 
acting fermion systems |L0]. In the phenomenological theory the coupling between electrons is 



given in terms of a few parameters which are determined by comparison with experiment [11 



Landau Fermi liquid theory, focusing on the collective excitations of the Fermi surface, is a 
type of bosonization. Zero sound modes, ordinary sound, and paramagnetic spin waves all 
have a bosonic character, despite their origin in terms of fermions. 

A microscopic interpretation of the Landau theory was provided in the late 1950 's and 
early 1960's via many body perturbation theory [|T2|, |X3 1 . Although the main idea behind 



Fermi liquid theory is deceptively simple, the concept of a quasiparticle, its realization in 
the microscopic theory is far from straightforward. Qualitatively when an electron interacts 
with all of the other electrons it polarizes its immediate vicinity due to the electron-electron 
repulsion, and this complex entity is the Landau quasiparticle. In this picture the interaction 
between electrons does not affect the quantum numbers of the electron, only its dynamics. 
For example the effective mass m* of the quasiparticle generally differs from the bare mass 
of the electron. In the microscopic approach, the existence of the quasiparticle can be 
ascertained by identifying a simple pole in the single-particle Green's function. At zero 
temperature this pole yields a ^-function peak in the single-particle spectral function at the 
Fermi surface. In the non-interacting case all of the weight is in this peak and there is 
no incoherent contribution to the spectral function. For free electrons the single-particle 
occupancy has a discontinuity of magnitude Zp = 1 at |k| = hp. When interactions are 
turned on, according to the Landau hypothesis, the discontinuity remains, but is reduced in 
magnitude to < Zp < 1, as shown in figure [I] with the remaining spectral weight appearing 
as incoherent background. This idea is made explicit by considering a power-series expansion 
of the self-energy of the quasiparticle about the low-energy limit: 

k) = a u + b e u + ic u 2 sgn(uj) + 0(uj 3 ) . (1.19) 

In the absence of singular interactions, Pauli blocking constrains the phase space for scatter- 
ing, as shown in figure 0. The scattering rate vanishes as uo 2 in the low-energy limit [14, [15 



and thus the £7(l)/f symmetry is recovered in this limit (see below). The mapping of the 
bare fermions Ck onto the Landau quasiparticles ipk may also be made explicit by writing: 

^ = Z^c k (1.20) 
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F 

Figure 1: Particle occupancy of the interacting Fermi liquid. The discontinuity at the Fermi 
surface has magnitude Z < 1; remaining spectral weight shows up as incoherent structure 
at higher energies away from the Fermi surface. 




(a) (b) 

Figure 2: A quasiparticle (a) of energy e can lower its energy via the formation of a particle- 
hole pair (b), but the phase-space for these processes diminishes as e 2 in the e — > limit. 
This is known as Pauli blocking. 
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where in terms of the self-energy, equation 1.19 , we have Z F X = 1 + a and m* = m x (1 



a)/(l + b). With this renormalization, the single-particle occupancy of the quasiparticles, 
(V'kV'k)) like that of free fermions, has a unit step at the Fermi surface. According to 



Luttinger's theorem ||16||, the discontinuity in the occupancy encloses the same volume in 
momentum space. In the special case of the isotropic Fermi gas with a spherical Fermi 
surface, the surface remains fixed at exactly the same position. 

The existence of Landau quasiparticles implies the existence of a Fermi surface, but the 
reverse does not necessarily hold. In one spatial dimension it is well known that the Landau 
picture breaks down, perturbation theory diverges, and the quasiparticle is destroyed by the 
interactions. In this case a different paradigm, the "Luttinger liquid," (the term having been 



introduced in 1981 by Haldane|FjJ) describes a wide range of interacting fermion systems. 



The foundation for this picture was provided in early work by Tomonaga[18[ who showed 
how to bosonize a one- dimensional interacting fermion system with a band cutoff. The foun- 
dations of one-dimensional bosonization were further developed by Luther and Peschel fT9[| , 



Coleman p0|) an d Mandelstam^TJ]. Mattis and Lieb|22| showed that a model of interacting 



fermions proposed by Luttinger p5[ (a model which allowed only forward scattering and in 
which the cutoff was taken to infinity) could be solved exactly. This solution demonstrated 
explicitly that the low energy degrees of freedom of the interacting one-dimensional electron 
gas are collective modes of charge density and spin density. Nonetheless a Fermi surface 
survives in the sense that the single-particle occupancy is non-analytic at exactly the same 
locations in momentum space as the free Fermi gas, 

nk ~ nk F — const sgn(/c — kp) \k — /c_f| 2< ' , (1-21) 

where the exponent ( depends on the interaction strength^. The Luttinger liquid fixed point 
has U(1)l ®U(1)r symmetry which reflects the conservation of the number of quasiparticles 
at each of the two Fermi points. Due to this enhanced symmetry, the Hamiltonian can 
be expressed in terms of left and right current operators, Jl{x) = : ip Jf L (x)il> '■ and 
Jr{x) = : ip R (x)ipn(x) :, and the bare two-body interaction can be replaced in the low- 
energy limit by the effective interaction 

H int = Jdx J L (x)J R (x) . (1.22) 

As the current operators obey a simple type of harmonic-oscillator algebra, the problem 
reduces to that of diagonalizing a Gaussian Hamiltonian. 

In spatial dimension D > 1, Haldane pointed out that the U(1)l <8> U(1)r symmetry 
enlarges, in the case of a Landau Fermi liquid, to the infinite U(l)'j? symmetry]!]. Although 
the bare interacting fermion Hamiltonian has only global £7(1) symmetry reflecting global 
conservation of charge, the Landau semiclassical energy functional 

e k = 4 + E/kk' 5n k , , (1.23) 

k' 

where Sn^ is the quasiparticle occupancy and k and k' denote momenta on the Fermi surface, 
has the same U(l)^ symmetry as the free system because the energy ey_ depends only on 
the occupancies n k , each of which separately is invariant under the U{1)^ transformation 
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equation 1.10 . Haldane proposed that a fully quantum mechanical description of the inter- 
acting fermion system might be obtained by bosonization, at least in cases where there is 
a well defined Fermi surface. In this viewpoint the formation of the Fermi surface is a zero 
temperature quantum critical phenomenon, an idea anticipated by Anderson in his book 
Basic Notions of Condensed Matter Physics ^ZQ 
exhibit the enlarged ^7(1)^° symmetry; 



Of course the bare Hamiltonian does not 
rather only the effective low-energy theory has this 
symmetry. Symmetry-breaking terms are irrelevant perturbations of the fixed point effective 
Hamiltonian. Indeed, for fermions with short-range interactions, this was shown explicitly 
in renormalization group (RG) calculations in a mathematical-physics style by Feldman and 
Trubowitz|25[ and by Benfatto and Gallavotti|[26||. Shankar's particularly clear exposition 



of the RG approach to interacting fermions is the one we draw upon most heavily in the 
following [p7j . 



A pioneering attempt to generalize the bosonization approach to spatial dimension greater 

Luther studied free fermions in general dimension by 



than one was made by Luther [28 



bosonizing the fermion fields at each point on the Fermi surface. The central idea was that 
the spectrum of the free fermion Hamiltonian can be reproduced by the particle hole excited 
states in the radial direction in momentum space, 



H = V F XXl k l ~ M V>kVka 
k 



1.24) 



and therefore it is possible to take the fermions to be one dimensional in each radial direc- 
tion k. In this construction, however, conservation of momentum makes it impossible for 
particle-hole pairs in different radial sectors to interact. This is where the subject remained 
until Haldane H showed that coupling between different points of the Fermi surface could 
be incorporated by adopting a different geometry: he coarse grained the degrees of free- 
dom near the Fermi surface into patches of finite extent rather than along radial rays. The 
approach was developed as a calculational tool by the authors and their collaborators who 

j. 133]. Independently, Castro-Neto 



emphasized the role of the current algebras [29, 30, 31 



and Fradkin used a coherent-state formalism which emphasized fluctuations in the geometry 
of the Fermi surface to derive similar results |]34], |35|. A program to include systematically 



corrections to the Gaussian bosonic fixed point has been initiated by Frohlich, Gotschmann, 
and Marchetti ||36|| and also by Kopietz and co-workers [37 , |38| . One of the goals of this ap- 



proach, which employs a Hubbard- Stratonovich transformation to replace the fermions with 
bosons, is to include non-linear terms in the fermion dispersion, as these terms are especially 
important in gauge theories of interacting fermions. See also work by Khveshchenko p>9" 



Some of these later developments are summarized in a review article by Metzner, Castellani 
and DiCastroffKJ] who discuss the close connection between multidimensional bosonization 
and the Ward identity approach pi] . Efl . 

It has been demonstrated how the main results of Fermi liquid theory can be recov- 
ered within this multidimensional bosonization approach^], ^ [Jl], [TT], [33], j3"7| . 
Furthermore, the stability of the Fermi liquid fixed point can be tested f30|, pT| • More impor- 
tant, bosonization provides a non-perturbative approach to the study of non-Fermi liquid 
fixed points that arise for example in the study of single layer quantum Hall fluids at even- 
denominator fillings v = 1/2, 1/4, . . ., and which may possibly be relevant in the cuprate 
superconductors, as we discuss below in Section [TO. In this instance, bosonization is a partic- 
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ularly appropriate tool, as within limitations described below it does not presuppose a Fermi 
liquid form for the quasiparticle propagator, and therefore non-trivial zero-temperature fixed 
points can be accessed. 

The outline of the rest of our review is as follows. In Section |2| we introduce the essen- 
tials of multidimensional bosonization. In Section |3] we employ the renormalization group 
to examine the stability of the Fermi liquid fixed point. The thermodynamics of interacting 
fermions is examined in Section [|. Collective modes are exhibited easily in the bosonized 
picture, more easily perhaps than in the traditional fermion approach, and an illustrative 
calculation of the spectrum of the charge mode is presented in Section The difficult prob- 
lem of calculating the fermion propagator is addressed in Section |6|. Two-particle properties 
such as the linear response to external density perturbations are examined in Section 
Although bosonization is limited to a description of low energy phenomena, high momen- 
tum processes are contained accurately. In particular, we show that the Kohn anomaly at 
wavevector |q| « 2k f is recovered within bosonization. In Section |8| we turn our attention to 
the non-Fermi liquid fixed points exemplified by composite fermion theories of the half-filled 
Landau level in the fractional quantum Hall effect. The gauge-invariant density response 
is calculated in Section |9|. We discuss nested Fermi surfaces and their possible relevance 
to cuprate superconductivity in Section IU. Finally we conclude in Section [TT] with some 



further thoughts about the usefulness of multidimensional bosonization. 



2 Basics of Bosonization 

In this section we formulate the problem of an interacting Fermi liquid using a bosonization 
scheme which is applicable in any spatial dimension. We start from the bare Hamiltonian 
of fermions interacting via two-body interactions. Initially, for simplicity, we consider a 
spherical Fermi surface. The bare Hamiltonian is given by 



H = Jd D x C t«(x) (^pj c Q (x) + IJ d D xd D y V(x - y) C t«(x) C ^(y) C/3 (y)c Q (x) . (2.1) 



First we make use of the renormalization group to integrate out the high-energy degrees of 
freedom [27 1 . The region of integration is depicted in figure ^|. The result is a low energy 



effective Hamiltonian describing quasiparticles with an effective mass m* interacting via 
two-body and higher-order interactions which, in general, include the long range Coulomb 
interaction as well as short range Fermi liquid interactions. Field operators ip k and ipk create 
and destroy these quasiparticles and are related to the bare fermion operators by 

^ k = Z k 1/2 c k (2.2) 

for momenta k which are restricted to a narrow shell of thickness A around the Fermi surface 
k F — A/ 2 < |k| < k F + A/2 depicted in figure [|. Note that an attempt is not being made to 
calculate the values of m*, Z^, or the other effective parameters from first principles; there 
is no generally reliable way to do this. Rather we focus on the remaining low-energy degrees 
of freedom and the effective theories which control them. 

As the next step, the Fermi surface is coarse grained by tiling it with squat boxes of 
height A perpendicular to the Fermi surface and dimension A ^ 1 along the surface. The 
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Figure 3: High-energy degrees of freedom, indicated by the hatched region, are integrated 
out yielding an effective theory for the remaining low-energy degrees of freedom near the 
Fermi surface (circle). These are to be bosonized. 
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Figure 4: The narrow shell of thickness A which contains the Fermi surface. The surface 
is tiled with squat boxes; the blowup shows one of the boxes at position S which contains 
a particle-hole pair of momentum q. We require |q| < A <C A <C kp for multidimensional 
bosonization to work (see text). 



boxes are small and have a squat aspect ratio in the relevant limits A <C A <C &p and 
A 2 /kp <C A. The reason for these limits, which are of crucial import, is discussed below. We 
introduce operators ip a (S; x) and ?/^ Q (S; x) which respectively annihilate and create fermion 
quasiparticles inside each patch 

</> Q (S; x) = e lks ' x £ 0(S; p)e i( -^> x ^ pa , (2.3) 
p 

where ks is the Fermi momentum located at the center of the squat box labeled by the 
coordinate S. Also, we introduce 6(S; p) which equals one if p lies inside the box and is zero 
otherwise. In terms of these operators the effective interaction is an instantaneous two-body 
interaction of the form 

V ta (S; x)^ t/3 (T; y)^(U, y)^ Q (V, x) . (2.4) 

This is the only marginal interaction as all terms with more operators or derivatives are 
irrelevant in the usual renormalization group sense |27]]; retardation may be neglected for 
the same reason. The bare Hamiltonian, equation |2.1|, is rotationally invariant, and barring 



spontaneous magnetization we can restrict our attention to rotationally invariant effective 
interactions. The interaction is heavily constrained by momentum conservation which de- 
mands that 

|ku + k v - k s - k T | < A < A . (2.5) 

In two dimensions there are only three possibilities: V = S and U = T (forward scattering); 
U = S, V = T (exchange scattering); and S = — T, U = —V (BCS Cooper pair scattering, 
where the (-) sign connotes the antipodal point). The three scattering processes are depicted 
in figure |5| and the corresponding Feynman diagrams appear in figure |6|. In three dimensions 
there is the additional complication of out-of-plane scattering |27[]. For the sake of simplicity 
we ignore these azimuthal processes here. 

It is clear that the interactions in the forward and exchange channels possess the large 
[7(1)£° symmetry: a U{1) phase rotation of both spin species of fermions in each squat box 
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Forward 



Exchange 



BCS 



Figure 5: The three marginally relevant scattering channels in D = 2 which are confined to 
the narrow shell of width A about the Fermi surface and which conserve momentum. 









(a) 



(b) 



(c) 



Figure 6: Feynman diagrams corresponding to the three marginally relevant scattering 
channels. The interactions between fermions (solid lines) are mediated by gauge fields (dotted 
line), (a) Forward scattering with amplitude / C (S, T) — | / S (S, T). (b) Exchange scattering 
with amplitude \ / S (S,T). (c) BCS Cooper pair scattering with amplitude Vbcs(S — T). 
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leaves the interaction in these two channels invariant as patches S and T each contain a 
creation and an annihilation operator. (There is no corresponding SU(2)f spin symmetry, 
as the exchange of spin between two separate patches on the Fermi surface breaks this 
symmetry.) The BCS interaction which couples quasiparticles in four separate patches breaks 
the U(l)f symmetry down to the ordinary global £7(1) symmetry which just reflects overall 
conservation of charge. We defer analysis of the BCS term until the next section, where it 
is shown to be marginally irrelevant or relevant depending on the sign of the coupling, and 
consider the remainder of the effective Hamiltonian: 

H[^M = /AE^(S;x)(^)i(S;x)} 

+ / d D xd D y ^{/ c (S,T;x-y)[^ ta (S;x)^(S;x)][^(T;y)^(T;y)] 

ST 

+ l/ s (S,T;x- y )[^(S;x)^^(S;x)].[^7( T;y )a>(T; y )]}. (2.6) 

The quasiparticle spectrum has been linearized about the Fermi surface, e(k) = ?/p(|k| — fop), 
and the quasiparticle velocity defines the effective mass through the relation v* F = kp/m*. 
Functions / C (S, T; x — y) and / S (S, T; x — y) parameterize the effective interactions in the 
forward and exchange channels. Upon making use of the Pauli matrix identity 

2 • o* = 28W - (2-7) 



it is straightforward to show that the forward scattering channel has amplitude f c — f s /4 
while the exchange channel has amplitude f s /2. Again no attempt is made here to calculate 
the short-range part of the parameters f c and f s as this is of comparable difficulty to a 
first-principles calculation of, say, the boiling point of water. Long-range density- density 
interactions such as the Coulomb repulsion, on the other hand, determine ||43|| the form 
of / c (S,T;x — y) in the limit of large separation, |x — y| — > oo. The remaining low- 
energy degrees of freedom within the shell are responsible for the screening of this long-range 
interaction, if screening occurs (see Section ||). 

To bosonize the effective Hamiltonian, equation ETB, we introduce the coarse grained 
currents 

J a (S; q) = £ f (S; k - q)0(S; k) {4_ q cA« " C^} (2.8) 

k 

here rik a = (^ka^ka) and the spin index a is not summed over, so all of the spin indices 
appear as subscripts. The index S labels a patch, with S = (9, 0) in three dimensions, on the 
Fermi surface at momentum kg. Recall that 8(S; k) = 1 if k lies inside a squat box centered 
on S with height A in the radial (energy) direction and area A D_1 along the Fermi surface; 
9(S; k) = otherwise. These two scales are made small in the sense that kp A ^> A > |q| 
by setting A = kp/N, and A = kp/N a where a < 1 and N — >• oo. We also require that 
a > 1/2 as otherwise the circular Fermi surface curves out of the squat box in the N — > oo 
limit. The relationship between these currents and the usual total density operator p(q) and 
total current operator J(q) familiar to condensed matter physicists is, for small-q, 

p(q) = E E 4.(8; q) (2.9) 
s <*=u 
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and 



J(q) = E E V HS) ^(s;q) 



(2.10) 



4-1 



where vp(S) = ks/m* is the fermion velocity at patch S. 

The coarse-grained current J Q (S; — q) creates a particle-hole pair of relative momentum 
q and spin a in patch S of the Fermi surface. As discussed in the Introduction, we expect 
Jq,(S; q) to have a bosonic character. Indeed, the currents satisfy the U(l) current algebra 



[J tt (S; q), MT; p)] = 5^ S^ 1 5° +pfi Q n s ■ q 



Here the quantity 



n = A 2J - 1 ( — 



2~ 



(2.11) 



(2.12) 



which equals the number of states in the squat box divided by the cutoff A, appears for 
the first of many times. Also rig is the unit normal to the Fermi surface at patch S. This 
result can be made plausible by the following argument. First, the current operators J a (S; q) 
and </g(T; p) clearly commute if a ^ (3. They also commute if S ^ T as the creation and 



annihilation operators are in separate boxes. The right hand side of equation [2.11| , if it is a c- 
number and not an operator, must be non-zero only when q = — p as otherwise momentum 
is not conserved. As the anomaly must change sign upon interchanging the two current 
operators, it must contain odd powers of q. Finally, as the current operators are scalars, a 
dot product of q with some other vector must appear. The only other vector is ns. 



The right hand side of equation 2.11 can be determined rigorously by direct computation 



on making use of the canonical anticommutation relations {iplf* , ippp} 
[J a (S;q),^(T;p)] 



t) a t) D 

°/3°q,p- 



£ 0(S; k - q - p)0(S; k) [0(S; k - q) - 0(S; k - p)] < +p>0 n ka 



+ £ 0(S; k - q - p)0(S; k) [0(S; k - q) - 0(S; k - p)] 

k 

^^^T l5 q+p,Q ^ n s • q + Error. 



(2.13) 



The first term on the right hand side of equation |2.13| , called the "quantum anomaly" or 
simply "anomaly," is of order XA d ~ 1 (L/2tt) d . The second term, the error, can be bounded 
by replacing the operator with a c-number: 



k— q— pa 



'1-5. 



q+p,oJ 



(2.14) 



and computing the volume of the geometrical intersection of the two 6 functions [0(S; k — 
q) — 0(S; k — p)]. Here a crucial assumption is introduced: the quasiparticle occupancy 
n ka = 1 deep below the Fermi surface on the lower face of the squat box, |k| = k F — A/2. 
Likewise, = on the upper face. Hence, matrix elements of the error term vanish on 
the upper and lower faces of the squat box, and a simple computation then shows that 



Error < max{|q|, |p|} x XA D - 2 {L/2n) D < X 2 A d -\L/2ti) 



D 



(2.15) 
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where L is the length of one side of the system. The error therefore is negligible compared 
to the anomaly as long as the limit A/ A — > (or equivalently N — > oo), is chosen. This 
establishes the key U(l) current algebra, equation |2.11| . The calculation of the error term 
makes it clear why the squat aspect ratio is so important: Scattering of quasiparticles between 
adjacent boxes is minimized by the construction^]. Unlike bosonization in one spatial 
dimension, which is accurate even when the cutoff A is an sizable fraction of the Fermi 
momentum, in higher spatial dimensions bosonization is more delicate. Because the Fermi 
surface curves, we require A <C kp, and to minimize interpatch scattering we require A A; 
hence A <C k F and the quasiparticle occupancy must change from 1 to over the small 
momentum range A. This point has not received sufficient attention in the literature. 
The charge currents also have an "Abelian" bosonic representation. If we set 

J a (S; x) = v 7 ^ n s • V0 a (S; x), (2.16) 

then we must demand that the bosonic fields a (S;x) satisfy the commutation relations: 



(S;x),^(T;y)] 



0: 



<W <*s,t e[ns ' ( x - y)]; l A s x (x - y)|A < 1 

|n s x (x-y)|A> 1 . 



(2.17) 



Here e(x) — 1 if x > and — 1 if x < 0. This can be accomplished by decomposing the boson 
fields into canonical boson annihilation and creation operators 



fS;x) 



i E 

q,n s q>0 



4(S;-q)e-^ x -a a (S;q)e^ 



where 



then 



;S;q),4(T;p 



v^s-q 



cD-l cD c 



a,/3 



J a (S; q) = V^|ns • q 



(2.18) 

(2.19) 
(2.20) 



t (S;q)0(ns-q) + at(S;-q)0(-n s .q) 

and equation ^ . 1 1| follows immediately. We have thus replaced the fermion bilinear appearing 
in equation [2.8| with a single boson operator. The boson operator acts on the quiescent Fermi 
liquid to produce a single particle-hole pair, which as noted above and in the Introduction, 
has bosonic character. 

We now introduce the key formula of the bosonization procedure which expresses the 
fermion quasiparticle field operator if> as an exponential ||44j| of the boson fields: 



^a(S;x) 




S;x) 6(S) . 



(2.21) 



Here V is the system volume, a = 1/A is the ultraviolet cutoff, and O(S) is an ordering 
operator introduced to maintain Fermi statistics in the angular direction along the Fermi 
surface 



O(S) = exp 



S-l 



ITT 



E E -4(T;q = o) 



T=la=T,J. 



(2.22) 
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Figure 7: The Fermi quasi-particle formed by repeated action of the operator <p which 
produces particle-hole pairs. One hole is left behind beneath the cutoff, so there are an odd 
number of fermions within the low-energy region specified by the ultraviolet cutoff A. 



where the patch labels T have been arranged consecutively. For example, the order can be 
chosen to start at point T = 1 on the Fermi surface, spiral outward in a clockwise fashion, 
and finally converge at the corresponding antipodep8| . The operator J Q (S; 0) appearing in 



equation |2.22| simply counts the number of quasiparticles in box T, so O(S) supplies the 



extra minus signs, much like the Jordan- Wigner transformation , which preserves anti- 
commuting statistics. In practice it is easy to identify these minus signs by hand and the 
explicit form of O(S) is not needed usually. Anticommutation rules are satisfied automati- 
cally in the direction perpendicular to the Fermi surface, as may be verified easily either by 
using the commutation relations, equation |2.17| , or by examination of the fermion propagator 



(see Appendix |C|). How can the fermion operator ip possibly be equated to the exponential 
of the bosonic operator 0? To understand this remarkable relationship it is important to 
appreciate the crucial role played by the cutoff A. As depicted in figure [7|, we can think of ^ 
as a sort of coherent state formed out of multiple particle-hole pairs||. A hole is left beneath 
the cutoff, and is therefore beyond the horizon of the low-energy degrees of freedom. 

It is important to test for closure using this bosonized representation of the fermion fields. 
We seek to verify that 

J a (S;x) = V hm{v4(S;x + en s )ip a (S; x) - (^(S;x+ en s )'0a(S;x)}) 

= v^ns- V0 Q (S;x) . (2.23) 

Here, in the first line, the current operator has been normal-ordered by point-splitting. That 
is, the expectation value of the fermion bilinear with fields at slightly different positions x 
and x + efis has been subtracted, to accord with the definition in momentum space, equation 
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1| To obtain the second line, first use the bosonization formula, equation [2.21|, to write 



V4(S;x + en s )V> a (S;x) = ^~ ex P 



Air 



—i- 



n 



> a (S;x + en s ) 



x exp 



47T 

«-^-0 Q (S;x) 



. (2.24) 



Then apply the following useful operator identity f46|j : 



e A e B 



: e 



A+B . (AB+ a2 + b2 ) 



(2.25) 



which holds if A and B are linear combinations of simple harmonic oscillator creation and 
annihilation operators. Here the notation : A : indicates normal ordering, with all annihila- 
tion operators on the right, and all creation operators on the left. This useful relationship, 
which we use frequently, can be derived using the more familiar Baker-Hausdorff relation: 



e A e B e -\ { A,B] 



(2.26) 



Now we can write 



V£(S;x + ens)V> a (S;x) = =p- lim : exp -i— ^ [0 a (S; x + en s ) - 0(S; x)] : 

4:71 

x exp— (0 Q (S;x + en s )0 a (S;x) -0 Q (S;x)) . 



(2.27) 



The equal-time correlation function of the <fi fields, derived in Appendix [A], is given by 

|n s x x|A < 1 



G^(S;x) 



^(S;x)0 a (S;O) 



(S;0)> 



n 2 j n ia . 

4-7T ns-x+ia' 



(2.28) 



— oo; 



Ins x x|A 3> 1 



Inserting this result into equation 2.27| yields: 



V£(S; x + en s )^(S; x) = iL : [1 - ie^n s • V0 Q (S, x) + 0(e 2 )] : 

Va si 



/a 



(2.29) 



The first, constant, term in this expansion about e = is subtracted off in the point- 
splitting convention of equation 2.23 and the continuum limit may be taken. The correct 
order of limits is lim e ^o lim a ->o because the separation between the fields e must always be 
much greater than the ultraviolet lattice cutoff a. As required by consistency J a (S;x) = 
^/Arrhs ■ V0 a (S; x). Similarly, by careful use of point-splitting it can be shown that the free 
part of the Hamiltonian is quadratic in the currents, 

Ho = v'f £ / d D x <jt°( S; x) (ii s • j) </>«(S; x) 



2ttv* f 



£ / ^[ns-V^Sjx)] 1 



(2.30) 
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quadratic because it has been assumed that the spectrum can be linearized near the Fermi 
surface. Non-linear corrections to the dispersion, when significant, lead to cubic and higher 
order terms in the boson fields which can only be treated perturbatively |T7[| . As a final 



check on the bosonization formula we compute the equal-time free fermion Green's function 



by again making use of equation 2.21 



G^S;x)HVaS;x)V4(T;0)> = 6^ e* k - x exp 



An 

— G^(S;x) 



Va 

1 tttd - ■ ■ ; n s x x A < 1 
0; otherwise 



(2.31) 



as expected. The low-energy fermion quasiparticle is related to the patch fermions by a sum 
over the patches 

^ a (x) = $>a(S;x) , (2.32) 
s 

so the usual Green's function Gp(x) can be obtained by summing over the patches, 

GV(x)=£'GV(S;x), (2.33) 
s 

where the prime indicates that the sum is only over patches for which |x x ns|A < 1. 
Physically this means that fermions in each patch, if)(S;x), propagate along nearly straight 
rays, reminiscent of semiclassical ballistic trajectories. 

The U(l)™ transformation rotates the phase of the fermions by a different amount 8s a 
in each patch: 

^ Q (S;x) ^e^ Q (S;x) . (2.34) 

The transformation leaves the current operators J a (S; x) invariant, as the phase factors ±#sa 
cancel out in the first line of equation |2.23| . Therefore Hamiltonians or actions which can be 



expressed solely in terms of the patch currents possess the large Z7(l)£° symmetry. In terms 
of the boson fields, the U(l)™ transformations amount to a translation of the fields by a 
constant: 

Q (S;x) ^0 a (S;x) + ^= 9 Sa . (2.35) 



47T 

Any Hamiltonian or action which conserves total charge is invariant under the global U(l) 
symmetry 6s a = The BCS channel breaks the U(l)^ symmetry down to this global 
symmetry, as the phases in the four separate patches (S, — S) and (T, — T) do not cancel out 
in general. The exchange channel, on the other hand, possesses a restricted U(l)^ symmetry: 
both spin species must rotate through the same angle, 9s a = #s- The restriction reflects the 
fact that while the exchange channel transports spin between two different patches of the 
Fermi surface, the total number of quasiparticles in each patch remains unchanged. Effective 
actions with only forward scattering possess the full U(l)f symmetry, however, and these 
are particularly amenable to bosonization, as shown below. 

Having shown how fermion quasiparticles are bosonized in each patch, we can now 
bosonize the effective Hamiltonian, equation |2.6|. If we choose z as the spin quantization 
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axis, introduce charge and spin currents 



J c (S;x) = J T (S;x) + J;(S;x) * 
J.(S;x) = J T (S;x)- ^(Sjx) ^ 

and the corresponding charge and spin bosons 

c (S;x) = T (S;x) + i (S;x) 

s (S;x) = T (S;x)-<MS;x) 1 J 

the effective Hamiltonian can be written in terms of the currents as 

HW\1>\ = Jd D x £ |^[J 2 (S;x) + J 2 (S;x)] 

+ JdPxdPy E { /c(S ' ^) 2 X " ^ MS; x) J C (T; y) + /a(S ' ^ " y) J S (S; x) J,(T; y) 
+ / s (S,T;x-y) [^}(S; x)^(S; x)^|(T; y)^ T (T; y)+ T<-| } (2.38) 
In terms of the fields it becomes 

H[<P c As] = Jd D x E^{[As-V0 c (S;x)] 2 + [n s -V^(S;x)] 2 } 

+ [ d D xd D y E |4vr/ c (S, T; x — y)[n s • V0 c (S; x)][n T • V0 c (T; y)] 



+ 4tt/ s (S, T; x - y) [(n s • V0 s (S; x)] [n T • V0 s (T; y)] 



\/47T 

I 2(-J cos^[0 s (S;x)-0 s (T;y)])}> (2.39) 



Although the spin and charge degrees of freedom separate in the Hamiltonian, equation 2.39 



it does not follow necessarily that the spin and charge excitations propagate at different 
velocities. 

It is remarkable that the charge part of the Hamiltonian is completely bilinear in the 
currents, or equivalently, the C fields, and therefore it is invariant under the full U(l)f 
symmetry and can be diagonalized exactly. In this instance the multidimensional bosoniza- 
tion scheme has transformed the four point fermion interaction into a problem of Gaussian 
bosons. Diagonalization is non-trivial, but as we show in subsequent sections, we can cal- 
culate the multi-point correlation functions exactly. The spin sector of the Hamiltonian is 
more complicated since the last, x-y, term in equation |2.39| , the coupling between the x- 
components and the y-components of the spin current, is a nonlinear cosine function of the 
spin boson fields. Consequently it breaks the symmetry and the equations of motion 

are nonlinear. Nevertheless the low-energy behavior of the theory can be analyzed by the 
RG method we discuss in the next section. 
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3 Stability of Fermi Liquids 



In this section we examine the stability of the Fermi liquid fixed point against BCS pairing 
in a weak-coupling RG calculation. It is well known that an attractive interaction leads to a 
superconducting instability, and here we show how this occurs in the bosonic representation. 
In contrast to mean-field theories, bosonization when combined with the renormalization 
group provides an unbiased means for analyzing instabilities, as it treats all of the interaction 
channels on an equal footing. 

For simplicity we consider the case of spinless fermions in two dimensions as the extension 
to three dimensions is straightforward. Consider also a circular Fermi surface to eliminate 
the possibility of nesting which could lead to charge or spin density wave instabilities in 
the forward and exchange channels. These possibilities are discussed later in Section 10. 
To further simplify the problem, initially we turn off the interactions in the forward and 
exchange channels, as we shall see that these channels do not alter the renormalization of 
the BCS channel at leading order, at least for the case of the circular Fermi surface. 

In the BCS channel a fermion in patch S on the Fermi surface is paired with the fermion 
in patch, — S, directly opposite, see figure |[ In cylindrical coordinates patches S and — S 
correspond to 6 and 6 + ix. To carry out the RG calculation we use the path integral 
framework, as we wish to rescale both space and (imaginary) time. The BCS action expressed 
in terms of the four Fermi fields is: 



S int [i/;*,i/;] 



(It 



d 2 x 

S,T 



Vbcs(S 



N*(0) 



V>*(-T; x)ip*(T; x)ip(S; x)ip(-S; x) (3.1) 



here iV*(0) = m* /(2ir) is the density of states at the Fermi energy and S and T are restricted 
to range over only half of the Fermi surface. As the fermions are spinless Vbcs must change 
sign under inversion so Vbcs{&) — —Vbcs(& + Also the Hamiltonian is Hermitian and 
therefore Vbcs(S — T) = Vbcs(T — S). The interaction can be written in bosonic repre- 
sentation using the bosonization formula, equation |2.21 . In doing so it is important to keep 
track of the order of the fermion operators. Formally the correct sign is set by the ordering 
operator O(S) but as noted in Section Q in practice it is easier to determine the sign by 
inspection of the fermion operators. With these remarks in mind the action is written as 



S, 



int I 



AA 

(2tt) : 



E 

ST 



dr 



2 Vbcs(S - T) 
d x — — exp 



N*(0) 



An, 



n 



ir(x)] 



(3.2) 



where we have introduced the combination 0s (%) 
usual 2 + 1 dimensional space-time coordinate x 



= (f)(S;x) 
= (x,r). 



+ ■ 



-S;x) and have defined the 



Because an exponential of the boson field appears in the action, equation |3.2| , it is most 
convenient to implement the RG transformation in real space. The RG transformation 
rescales only the component of length parallel to the Fermi surface normal, x» = ns ■ x. The 
perpendicular component, x± = ris x x, remains invariant so that the real space volume 
element d 2 x = dx»dx± becomes (sdx\\)dx± and dr —>■ sdr. The quadratic part of the action 
is invariant under this transformation provided that the boson field is invariant, <fi — > <p. 
Then it follows that the BCS interaction is marginal as expected. The usual renormalization 
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group procedure is to integrate out the fast parts of the fields 0(S;x). Here "fast" means 
high-energy. To be precise modes with momenta 



A 



A 



^<|ns-q|<2 



(3.3) 



are eliminated which in real space translates into integrating out the short distance pieces 
of the fields with 2a < |ng • x| < 2sa, where we recall that a = A -1 is the ultraviolet cutoff. 
To complete the transformation xn — > sxn, and r — > sr. Integration over the fast modes is 
accomplished as usual within the functional integral 



-S[4>(S;x)] 



n 



J2a<|nx-x|<2as 



L>0(T;x) e 



-So [01— Si, 



(3.4) 



As only the fast parts of the fields are integrated out it is convenient to split the <p fields 
into slow (f)' and fast h fields by writing <fi — $ + h. Then we can write the right hand side 
of equation |T4] as exp(— So [</>]) x (exp(— S int [4>' + h])) where () means an average over the 
fast h(T; x) fields only. The renormalized interaction is found perturbatively by expanding 
in powers of Vbcs an d at leading order 



(Sint [ft]) 



AA 
2^ 



V BCS (S - T) 



X- 



J S,T 

Ttt 



N*(Q) 



exp < i- 



n 



x 



(j)' T (x)} 



X 



exp i—^- [h s (x) - h T (x)] 



(3.5) 



To evaluate the restricted average () we make use of the propagator of the h fields (the 
derivation of the propagator is given in Appendix 0): 



(h(S; x)h(S; 0)) 



■In 



0; 



In 



Xn +iVpT+(2isa)sgn(T) 
xn+iVpT+(2ia)sgn(T) 



In particular 



Now use the fact that, for a Gaussian measure, 



(e ih ) 



e -W) 



A|x_l| < 1 
otherwise. 



(3.6) 



(3.7) 



(3.f 



and equation 3.5 becomes 



(SintW}) 



AX/s 
2rr 



E 

S,T 



dr / d 2 x 



Vbcs( 



N*(0) 



exp 



(x)-0' T (x)] , (3.9) 



where the prime on the integral indicates that the region of spatial integration covers only 
small distances. Then, upon rescaling space and time, we recover Sj nt [</>]. 
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The first nontrivial contribution to the renormalized interaction is at second order in the 
expansion in powers of the BCS interaction, Vbcs- This term is evaluated by making use 



of the variant of the Baker-Hausdorff formula, equation |2.25| , and the propagator of the h 
fields: 



AA ^ 



(2vr) 2 s 



J2 [dr f d 2 x VBCS ^ S ^j^p 5 ^ — ex P { i^Wsix) - M%)\ 



J S,T 




X Ai-oo ° J 2a m [v* F a + 2a sgn(a)] 2 + y\ [ ' ' 

where y» = fiu ■ y is the spatial integration variable and the factor of 1/A in the second line 
is the result of integrating over the perpendicular coordinate, y±. As s — > oo the highest 
power of s (in this case ln(s)) gives the leading contribution to the renormalization. The 
time integral over a may be evaluated by contour integration. The structure is such that of 
the two poles, one lies in the upper half plane and one in the lower, which is a direct result 
of the coupling in the BCS channel between patches U and — U. The remaining spatial 
integral is proportional to ln(s) and we find that the flow of the coupling function in the 
BCS channel is given by: 



dV BCS (S - T) A 



Vbcs(S - U)V BGS {U - T) . (3.11) 



d\n(s) (2Tr) 2 k F ^ 
For the circular Fermi surface considered here 

and the flow equation can be diagonalized by a Fourier decomposition of the BCS interaction 
into angular momentum channels: 

1-2-k At) 

^ e m V BCS {6) . (3.13) 

JO 2,71 

The flow equations decouple channel by channel, and are simplified further as Vi is dimen- 
sionless: 

= -Yl . ( 3.i4) 

dln(s) 4tt 1 ' 

This result, which is identical to that derived by Shankar in the fermion basis |27[], tells us 
that the couplings Vg, which are marginal at tree level, are in fact marginally relevant if 
negative, and marginally irrelevant if positive. The first case leads to the BCS instability, in 
the second case the Fermi liquid fixed point is stable, as indicated by figure |[ 

We should ask now, what would happen were we to include at the outset the other 
marginal interactions, for example the coupling between the charge- currents: 

•S^ = t72 £ I dT d2x d 'y ^( S ' T; x - y) J C (S; x, r) J C (T; y, r) . (3.15) 

V ST 
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Fermi Liquid Fixed Point 




BCS Instability 

<<(((< ) | < <<( ( — ^ v 

/ 

Figure 8: RG flows decouple, channel by channel. The Fermi liquid fixed point is at Vg = 0. 

Since the interactions are quadratic in the bose fields, they parameterize infinitely many dif- 
ferent Gaussian fixed points which do not flow under the RG transformations. Incorporated 
into the RG calculation they contribute only subleading corrections to scaling and do not 
influence the leading RG flows of the BCS interaction. For example, although bosons in 
different patches are correlated via the interaction, this correlation is only a fraction, A/fcp, 
of the correlations within the same patch, as we prove below in Section |6|. The leading RG 



flows, equation |3.14| , are unchanged. 

It is not difficult to see that the system develops an energy gap if the BCS interaction is 
attractive in any channel. The RG flows in attractive channels scale to large negative values 
as s — > oo, and the size of the gap may be estimated by calculating the length scale at which 
the renormalized coupling becomes of order one. Alternatively, a mean-field approximation 
shows how the gap develops. Suppose for simplicity that Vbcs(@) — Ve=iCos(9), conducive 
to p-wave ordering. As in the classic BCS calculation we may decouple the Cooper pairs via 
a Hubbard- Stratonovich transformation upon introducing a field ^(S) = A(S)e l *. Then the 
bosonized Hamiltonian for the pair (S, — S) is given by 

ff[0s] = /^{^T ([(ns-VW.(S;x)] 2 +[(n s .V)0(-S;x)] : 



At the mean field level, A(S) can be assumed to be constant with no fluctuations, and the 
problem is equivalent to a sine-Gordon model in each direction lis- The continuous U(l)f 
gauge symmetry of the action under translations of the boson field by an additive constant, 
0s — > 0s + ~^(^s + #-s), breaks spontaneously giving rise to the usual superconducting 



energy gap. 

There is in fact the possibility of a BCS instability even when the interactions are repul- 
sive: the Kohn-Luttinger effect |47|]. This possibility only becomes apparent beyond leading 
order. The bare couplings, Vg, due to a short-range repulsive interaction are all positive but 
tend rapidly to zero at large £. Fermi liquid interactions as discussed above generate irrele- 
vant contributions in the BCS channel down by a positive power of A/fcp, nevertheless this 
is sufficient to make some of the Ve slightly negative, at least at sufficiently large £, leading 
to an instability. Because of the small magnitude of these coefficients, the effect can only 
be important at extremely low temperatures and the essential physics is controlled by the 
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Fermi liquid fixed point. However, nesting instabilities can enlarge dramatically this effect, 



and we examine this situation in Section 10 



Finally, it might be thought that in general an instability must occur in the spin sector 
since the interaction is a cosine function of the spin boson. Indeed this does occur in one 
spatial dimension: the coupling between spin currents at the left and right Fermi points flows 
under the RG transformation [E8, |9L [50, l6|. If this were to be the case in D > 1 it would 



spell disaster for the multidimensional bosonization approach to Landau Fermi liquids since 
the Landau parameters in the spin sector should be constants parameterizing the Landau 
fixed point. However, on closer inspection it is clear that the spin couplings are indeed 
fixed. Repeating the weak-coupling RG calculations to second order for the x-y part of the 
interaction in equation 2.39| gives 
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dr 
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#(T;s)] 
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+ iv* F a) 2 



(3.17) 



In this case, unlike the case of the BCS interaction, the correction vanishes since the integral 
over (y\\ , a) does not encompass a simple pole, and the spin interaction is strictly marginal, 
at least within this second-order perturbative calculation. A deeper, and non-perturbative, 
reason for the exact marginal nature of the spin current coupling in D > 1 is given below in 
Section [|. 



4 Thermodynamics of Fermi Liquids 

In this section we use the bosonization formalism to calculate the specific heat of the inter- 
acting Landau Fermi liquid in three dimensions. We show that a well-known nonanalyticity 
in the specific heat, a term proportional to T 3 lnT, is reproduced by the method. The ex- 
istence of such a term is consistent with careful measurements pTf of the specific heat of 
helium-3. It has also been seen in heavy fermion systems. For the sake of simplicity we turn 
off the spin-spin interaction and eliminate the spin index. Then the bosonic Hamiltonian of 
the interacting spinless fermion system is given by 

H = \ EE V(S, T; q) J(S; -q) J(T; q) (4.1) 
ST q 

where 

nS,T; q ) = | 5s V + ^^. (4.2) 

As the nonanalytic behavior of the specific heat derives from small momentum processes it 
is sufficient to take the Fermi surface to be locally flat. If we again let q\\ = fig • q denote the 
component of the momentum q parallel to the surface normal or Fermi momentum in patch 
S, the charge currents J(S; q) can be represented in terms of canonical boson creation and 
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annihilation operators as follows 



J(S;q) 



-Slq\\ a f (S; -q), q\\ < 



(4.3) 



Slq\\ a(S; q) 



> 



With the replacement of the current operators with the canonical bosons the Hamiltonian 
becomes 

H = E E ^( S > T; q) ng,| a t(S; q)o(T; q) (4.4) 

S,T q,9||>0 

which we can diagonalize directly to find the spectrum. The specific heat is then computed 



using the standard formula for bosons, equation 1.14 



First we review the case of noninteracting fermions, /(S,T;q) = 0, discussed in the 
Introduction. The eigenenergies of the bosonized Hamiltonian in this case are simply 



e(S; q) = v* F n s • q = v* F q\\ 



(4.5) 



which agrees with equation ETfil the energy of a particle-hole pair. The sum over the patch 
index S and the components of q perpendicular to the Fermi wavevector counts the number 
of states at the Fermi surface which we call A. The sum over the components of q parallel to 
the Fermi wavevector can be converted to an integral. Then assuming that the temperature 
is sufficiently low such that the thermally excited particle-hole pairs lie within the cutoff, 
k B T <C v F X, we obtain 
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Ak B T 2 



L 
2^ 



vfq\ 
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(4.6) 



This is the correct result for spinless fermions in D = 3 which, of course, should be multiplied 
by a factor of 2 to account for spin. It is remarkable that the boson formula, equation |1.14| , 
yields the full specific heat. Multidimensional bosonization reproduces the entire fermion 
Hilbert space. 

Now we follow Pethick and Carneiro[|5*2| and focus on particle- hole pairs separated only 
by a small momentum w = k$ — k T (with |w| <C hp) since a consideration of these processes 
is sufficient to demonstrate the existence of nonanalytic contributions to the specific heat. In 
this limit we may continue to approximate the Fermi surface as locally flat as shown in figure 
H this simplifies considerably the problem of diagonalizing the quadratic boson Hamiltonian 
(see Section ||). Now the current-current coupling /(S, T; q) may be expanded in powers of 
the dimensionless, and rotationally- invariant, expansion parameter q 2 /(w 2 + q 2 ): 



/(S,T;q) = a + b 



q 2 + w 2 



a + b 



n 



+ 



+ 



(4.7) 
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Figure 9: The Fermi surface treated as locally flat, with nearby particle-hole pairs of 
momentum q and — q coupled together. 



In the second line we have dropped the perpendicular q^ = fis x q terms in the numerator 
since they do not contribute to the non-analyticity. The same approximation can be made 
in the denominator, as the squat aspect-ratio of the patches enforces |w| ^> |q| in the 
iV — > oo limit. The q 2 term in the denominator must be retained to eliminate the artificial 
divergence which would otherwise occur in the w 2 — > limit. The expansion is controlled in 
low-temperature limit which ensures that |gn|, the particle-hole momentum perpendicular to 
the Fermi surface, is small. The current-current interaction coefficient b, which therefore is 
irrelevant in the RG sense, gives subleading corrections to the specific heat. In actual fact, 
the interaction considered here differs from that studied by Pethick and Carneiro in that it 
couples particle-hole pairs at point S and T whereas the Pethick-Carneiro interaction couples 
the quasiparticle occupancies n p and n u+p . To be precise, the Pethick-Carneiro interaction 
has the form b' E u .p(u'P) 2 ^p Sn u+p . This interaction cannot be expressed directly in terms 
of the currents since it involves the products of quasiparticle occupancies above and below the 
Fermi surface, whereas the current operator evaluated at zero momentum J(S; 0) averages 
the occupancy operator over the interior of a squat box. Therefore a direct connection with 
the earlier calculation of Pethick and Carneiro cannot be made. 

To proceed, we diagonalize the Hamiltonian with the aid of a Fourier transform from 
patch index space S to X space. Let 



t (S;q) = -^L= f d 2 X e" iX - s <J(X,q) 
■ kirk 2 , J 



(4- 



where / d 2 X 1 = Arckp / A 2 is equal to the number of patches which cover the Fermi surface. 
Then 

H = L 3 j d 2 X / A/2 ^ d ^Lf Xd ± e(X; q)at(X; q)a(X; q). (4.9) 

J J-A/2 Z7T J-A/2 2n Jo 2ir 

Using equation \i.'2\ and the representation of the current given in equation |4.3| we obtain the 
eigenenergies 

bA 2 q 3 e -«-w/A 



e(X; q) = v* F q\\ + 



E 



(4.10) 



The sum can be converted to a Riemann integral by the substitution A 2 E w ~" > I d 2 w and 



we find 



e(X; q) = v* F q\\ 



bqf. 



(27T) 



-vrln(g |X| 2 ) 



(4.11) 
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In this equation we have discarded those terms proportional to b that make analytic contri- 
butions to the specific heat and have retained only the logarithmic correction. We treat this 
term as a perturbation and calculate the specific heat to 0(b). The change in the specific 
heat 5Cv due to this perturbation is 
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(4.12) 



In the second line of this equation we have retained only the term containing T 3 In T tem- 
perature dependence. The integral in the second line equals 7r 4 /30 so the final result is 
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15k 2 F 



-T 3 InT 



(4.13) 



Not surprisingly this result has the same form as that found by Pethick and Carneiro[^] as 
dimensional analysis makes this likely. A direct comparison of the coefficient is meaningless, 
however, as the interactions are not the same. Nevertheless, the same point is made: sublead- 
ing corrections to the linear specific heat have many different origins. Irrelevant interactions 



such as equation [4.7| or the Pethick- Carneiro form both make nonanalytic contributions to 
the specific heat, as do other more complicated irrelevant terms |53|. Repeating the calcu- 
lation in D = 2 we find that the interaction equation [4.7| leads to a subleading correction 
5Cy ~ bT 2 , also in agreement with fermionic calculations p3 . 



5 Collective Modes 

The curvature of the Fermi surface did not play an important role in the determination of 
the specific heat. In fact we took the Fermi surface to be flat and consequently the Hamil- 
tonian could be rewritten as the sum of products of a single creation operator and a single 
annihilation operator. Collective excitations of the Fermi surface, on the other hand, are 
possible because of the curvature. The Fermi surface itself may be viewed as dynamical, 
sloshing back and forth relative to its equilibrium position in momentum space. It is im- 
portant to reproduce the collective mode spectrum within the bosonic representation. For 
a curved Fermi surface the Hamiltonian contains the products of two annihilation operators 
and the products of two creation operators in addition to the bilinear combination of anni- 
hilation and creation operators found for a flat Fermi surface. Therefore a straightforward 



diagonalization of the bosonic Hamiltonian is no longer possible, see however |34|, pq| . 

In order to discuss the spectrum of collective modes it is necessary to elaborate further 
on the Hamiltonian. In the absence of an attractive BCS interaction the excitations of the 
interacting Fermi liquid are in the particle-hole channel and carry either charge or spin. To 
exhibit this factorization into charge and spin sectors explicitly the bosonized Hamiltonian 
can be written as H = H c + H s . The Hamiltonian for the charge sector in D dimensions, as 
before, is bilinear in the charge current operators J(S; q) = J2 a Ja(S;q). We drop the "c" 
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subscript on the charge currents as no ambiguity results. 

H c = lj2J2 Vc(S, T; q) J(S; -q) J(T; q) (5.1) 
st q 

The prefactor of 1/2 compensates for the double-counting q <-> — q in the sum over the 
momentum. For the present we consider only the case of short-range Fermi liquid interactions 
/ C (S,T): 

V C (S, T; q) = ifl-^gi 1 + ~/ c (S, T) . (5.2) 

The charge currents satisfy the equal-time U(l) current algebra, equation 2.11| , but with 
twice the anomaly as the spin index is summed over. The factor of \ in the first term of 
equation FT2] compensates for this doubling. 



For our analysis of the spin collective modes it is convenient to introduce spin currents 
which are manifestly SU(2) covariant instead of continuing to employ the spin quantization 
axis as in Section [|. These so-called "non-Abelian" currents [|55| are defined to be: 

M q) = ^E *(S; k - q)0(S; k) ViV&kk (5-3) 

where again a are Pauli matrices. Note that, unlike the case of the charge currents, there 
is no need to subtract the vacuum expectation value as the quiescent Fermi liquid is non- 
magnetic. Like the charge currents, spin currents in different patches on the Fermi surface 
commute. They also commute with the charge currents. Spin currents in the same box 
satisfy the more complicated SU(2) current algebra: 

J a (S; q), J b (T; p)] = i Q n s ■ q c^ 1 5 ab 5° +P ,o + < T V bc J C (S; q + p) ; (5.4) 

here (a, b, c) = (x, y, z) label the three components of the spin, and e abc is the completely 
antisymmetric tensor. The coefficient of the first term on the right hand side of the cur- 



rent algebra, equation |5.4j , is known as the "central charge," k. In gapless one-dimensional 
antiferromagnets[5*6"] with half-odd-integer spin- S, generically k — 1 except at special inte- 
grable points where k = 2S. Here, by contrast, k = (AL/2tt) d ~ 1 ^> 1. The effective spin 
is large because the coarse-graining procedure sums over a multitude of momentum points 
inside each squat box; consequently the spin currents are semiclassical in nature. 

The spin Hamiltonian is quadratic in the non-Abelian spin currents (unlike the Abelian 
case, there are no cosines) and Fermi liquid spin-spin interactions f s which appear in the 
Hamiltonian couple together spin currents in different patches: 

= ^( s > T ; q) J ~( s ; -<*) ■ W q) ( 5 - 5 ) 

ST q 

where 

V S (S, T; q) = I^O-^gi 1 + i/.(S, T) . (5.6) 

Note the factor of 2/3 appearing in the construction of the free Hamiltonian, which we 
discuss below. 
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Instead of diagonalizing H directly, we find instead the equations of motion for the charge 
and spin currents which we shall show are equivalent to the well-known collective-mode 
equations. In the charge sector, 

iJ^(S;q) = [J(S;q),F c ] 

2Q 

= v* F h s -q J(S; q) + n s -q-rE /c(S, T) J(T; q) . (5.7) 

v T 

The first term on the right hand side of equation |5.7| is due to the free dispersion relation 
for the particle- hole pairs of momentum q in the patch S. The second term couples currents 
in different patches. Note that 

£E - /(g^l-M 

P /•! /-27T 



d(cos#) / # (£> = 3) ; (5i 



(2tt) 3 7-i v ' Jo 

so the second term reduces to the usual integral over the Fermi surface in the continuum 
limit of a large number of patches, N — > oo. The equations of motion for the spin currents 
contain, in addition to these two terms, a third term which makes the spins precess in 
the local magnetic field, a phenomenon which is especially important when the system is 
magnetically polarized. 

^J a (S;q) = 4ns-qJ fl (S;q) + |n s .q^/ s (S,T)J a (T;q) 

- ^ c £j 6 (S;k)£/ s (S,T)J c (T;q-k) . (5.9) 

V k T 

Note that the factor of 2/3 multiplying the Fermi velocity v* F in the free part of equa- 



tion |5l] does not appear in equation |5T9| . The origin of the factor of 2/3 is easy to un- 
derstand in the Sugawara construction of the free fermion Hamiltonian out of current 
bilinears|57|, |58fl . It reflects the SU(2) invariance of the spin currents which permits the 
replacement J(S; q) • J(S; — q) — > 3J 2 (S; q) J Z (S; — q), at least for the purpose of computing 
the spectrum. The derivation given here, on the other hand, does not rely on this argument 
as spin rotational invariance is respected explicitly. Rather the factor of 2/3 in equation |5.6| 
cancels contributions to the equations of motion that arise from both the 5 ab and ie abc terms 
in the spin current algebra. 

Now we are in a position to understand the deeper reason for the absence of renormaliza- 
tion of the spin Landau parameters alluded to in Section fj. Unlike the D = 1 case, the spin 
current operators defined by equation |5.3| involve a sum over a large number of momenta 



inside the squat box. This is apparent in the non-Abelian algebra equation |57| as the diag- 
onal anomaly (the first term proportional to S ab ) is of order 0| fig • q| ^> 1. In the absence of 
macroscopic spin polarization the anomaly overwhelms the off-diagonal, ie abc , term. Physi- 
cally this means that, in zero applied magnetic field, the spin equation is identical in form to 



the charge equation, and spin waves propagate freely as expected|TT|]. In the opposite limit 
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of a large external magnetic field, the last term in equation |5l| dominates. Spins precess 
locally in this limit. Returning now to the question of the spin Landau parameters, we see 
that they remain fixed, as only the non-linear ie abc term can drive renormalization flows and 
it can be neglected. This proof, which is non-perturbative, and therefore independent of the 
size of the spin Landau parameters, is consistent with the absence of renormalization found 
in the weak-coupling calculation presented in Section |3|. 

As it stands these operator equations are exact, at least in the iV — > oo limit in which 



the current algebras, equations |2.11| and |5.4| , become exact. Solving the charge collective 
mode equation |5.7| is straightforward as it is linear. The equation for the spin currents is 
more difficult to solve due to the nonlinear term. Here for completeness we solve the charge 
equation for a spherical Fermi surface in D = 3 with one Landau parameter, / C (S, T) = / . 
As usual we search for a normal mode solution of the form: 



J(S;q,i)=e-^u(S;q) 



(5.10) 



also set q = qz for definiteness. Then equation 5.7 becomes: 



(vtqcosd -u)u(9,<p;q) + v* F qcos9^ E - C d(cos0') dfiuW, </>'; q) = (5.11) 

Att 6 v*p 7-i Jo 



and the solution then is given by 



u(9,<f>;q) 



cos 9 
x — cos 9 



(5.12) 



where we have introduced the dimensionless quantity x = ^rj^| • Equation |5. 11| determines the 



allowed values of x. Substituting equation |5.12| into equation 5.11 we find that propagating 
modes exist for x > 1 with x given implicitly by: 



k F m*f f xln X + l_ 2 



2tt 2 



x — 1 



^ln^-2 
2 V x-l 



(5.13) 



with the introduction of the conventional dimensionless form of the Landau parameters 

F£ = iV*(0)i£ (5.14) 

and recalling that in D = 3 



^ (0) = 2 f-A^dk 8m _ eF) m*k F 



o (2vr) 3 



7T" 



(5.15) 



Analytic solutions of equation |5.13l exist HI 51 in the extreme limits < Fq 1 and Fq ^> 1: 

l + 2e~ 2 / F °, 0<F <1 



x 



Fq/3, 



(5.16) 



F > 1 



31 



It is remarkable that the charge sector can be solved exactly without taking the semiclassical 
limit of a macroscopically occupied collective mode. We show in the next Section that the 
position of the poles of the boson correlation function agree with the dispersion of the 
collective mode found here, as must be the case. 

To complete the discussion, we show how to construct spin collective modes in the semi- 
classical limit. We seek to replace the quantum spin current operator J a with its semiclassical 
expectation value: 

(J a (S;q,t)) = S a (S;q,t) (5.17) 

— * — # 

Note that the collective mode amplitudes are real valued in x-space as [J(S; q)]> = J(S; — q). 
The replacement of the current operators by c-numbers can be accomplished formally by 
introducing a coherent state basis which spans the space of geometric distortions of the 



Fermi surf ace |34], Spin collective mode coherent states are generated by exponentials of 



the spin current operator 



MS}) = exp E 2 S(T; q) • J(T; -q) j> |0> (5.18) 

I T q " n T - q J 

where |0) represents the quiescent Fermi liquid. A simple computation using the spin current 



algebra equation [5]4| then shows that 

ms}\nT;p)ms}) 

(*[§\M$\) 



S a (T;p) (5.19) 



consistent with our definition, equation |5.17| . In the equation for the spin collective modes it 



is also necessary to decouple the expectation value of the product of two spin currents into 
the product of expectation values S(T; k) x S(T; q — k). The decoupling is exact only in 
the semiclassical limit of a macroscopically occupied spin mode. 

In summary, multidimensional bosonization captures all of the well-known physics of 
charge and spin collective modes. The derivation of the equations of motion is straightforward 
and relies only on the existence of the Abelian and non-Abelian current-current commutation 
relations and the quadratic Hamiltonian which describes the £7(1)£° fixed point. 



6 Quasiparticle Properties 

In this section we use the bosonization scheme to determine the fermion propagator of a 
system of spinless fermions interacting via two body forces. The question we address is: 
What is the nature of the ground state and how does the ground state depend on the spatial 
dimension D and the range of the interactions? 



Making use of the bosonization formula, equation 2.21, the fermion propagator for a 



single patch on the Fermi surface can be expressed in terms of the boson propagator 
G F (S;x,r) = -<^(S; x, r)^(S; 0, 0)) 

= "(^jTi^ - (exp[i-^-0(S;x,r)]exp[-i-^-0(S;O,O)]) (6.1) 
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which can be rewritten by means of the Baker-Hausdorff formula, equation |2.26| as 



G F (S;x, t) 



h D ~ 1 e* ks ' x 



exp 



— (0(S;x,r)0(S;O,O)-0 2 (S;O,O)) 



(6.2) 



(2tt) d a 

It is convenient to express the fields, <p i n terms of the canonical boson fields, a and a*, 
equation |2.18| . Then the correlation function G</,(S; x, r) can be written as 



G (S; x, r) = (0(S; x, r) 0(S; 0, O)-0 2 (S; 0, 0)> = — £ / — [ e ^—)-l] 



(a(S;q)a f (S;q)) 



2tt 



n s ■ q 



(6-3) 

which determines the fermion Green's function once the Fourier transform of the boson 
Green's function 

G B (S; x, r) = -< a (S; x, r)a\S; 0, 0)) (6.4) 

is known. This procedure although straightforward in principle can be difficult to carry out 
in practice. As we show below, the Fourier transform of the boson Green's function can 
be obtained exactly for any interaction and for any spatial dimension, but to determine the 
fermion Green's function this result must be back transformed, equation |6l| and exponenti- 
ated, equation |6.2| . Finally, to extract the nature of the fermion propagator, this space-time 
result must be Fourier transformed back into momentum and frequency space. 

To determine the boson propagator we introduce source fields £(S; q) and £*(S; q) where 
q = (q, uS) is shorthand notation for the D+l dimensional vector and construct the generating 
functional: 



Z[Z,C]=jDa(S;q)Da*(S;q)exp -S[a*, a] - £ £ J {C(S; q)a(S; q) + £(S; q)a*(S 

(6.5) 



Q S,n s q>0 

Here the action is given in terms of the Hamiltonian as: 



S[a*,a]= J dr jy^ J d D x a*(S; x)-^-a(S; x) + H[a*, . 



(6.6) 



Then the boson Green's function may be found by differentiating with respect to the source 
fields: 

G B (S,T;g) = -(a(S;gK(T;g)>: ' >_ ' lnZ 



6tiT\q)6?(S;q) ' (6 ' ?) 

To illustrate this procedure consider a system of spinless fermions in two dimensions inter- 
acting via short range Fermi liquid interactions. The action is 



S[a*, a] = So [a*, a] + Sfl[o>*, a] 



(6.1 



where the free fermion contribution to the action, written in terms of the canonical boson 
fields, is: 



S [a*,a] = -J2 E f J~( iu ~ v *F h s ■ q)o*(S;g)a(S;g) 



(6.9) 



Q n s q>0 
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and the term due to Fermi liquid interactions can be expressed as: 

^L[o*,o] = ^/^EE/(S,T)J(S;-g)J(T; g ) (6.10) 

where we recall that 

J(S; q) = Jnh s • q [a(S; q)0(S; q) + a*(S; - g )0(S; -q)] . (6.11) 



In equation [6J| and in what follows we abbreviate E q Es,n s q>o by simply Z) q Z)n s .q>o- 
For free fermions the generating function is 



Z[£,C] = [DaDa*explj2 E / 7T a *( S 'i<l){ iuj ~ v *F&s • q)o(S; g) 



q n s q>0 

r(S;g)o(S;g)-e(S;g)a*(S;g) 



(6.12) 



The action is quadratic in the bose fields which can be integrated out to give 

^,C]=exp-£ £ f^e(S;q) T . ^ r£(S;g) (6.13) 

and therefore 



G B (S,T;q) = -{a(S-q)a*(T ]q )) 



D-l 



XL/— 
°S,T 



iu - Vpiis ■ q 



(6.14) 



On taking the modified Fourier transform, equation pT3| , the patch-diagonal piece of the 
correlation function is given by: 



G°(S;x,r) = -^ln 



n s ■ x + iv* F r 



id 



(6.15) 



Exponentiating this equation, the free fermion Green's function is found to be given by 



G° F (S;x,r) 



iA e iks ' x 
(27r) 2 h s ■ x + ivpr 



(6.16) 



Additional details of the calculation are presented in Appendix |C[ 

The analysis can be extended now to include the interactions, equation |6.10| . First Fourier 
decompose the coupling /(S,T) and introduce Landau parameters f n : 



/(S,T)=£/nCos[n0(ks,k T ) 

n 

and recall the addition theorem in D = 2 



(6.17) 



cos[n6>(k s , k T )] = (— l) n cos[n6>(k s , q)] cos[n#(k T , q)] + sin[n0(k s , q)] sin[n6>(k T , q)] 



(6.18) 
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Here 0(ks,k T ) is the angle between k s and k T . To avoid the complicated Bogoliubov 
transformation needed to diagonalize the Hamiltonian, we introduce instead auxiliary gauge 
fields A e n and A^, two fields for each n, and perform a Hubbard- Stratonovich transformation 
to decouple the interaction^]. We call these auxiliary fields "gauge" fields since in the case 
of the long-range Coulomb interaction which we consider below they correspond to the scalar 
component of the electromagnetic potential, A , which mediates the interaction (see figure 
Bf). Written in terms of the gauge fields the short-range interaction has the form 



S 



FL 



duo 



2V J 2tt ? 



-1)" Ai(-q)±-Ai(q) + AH-q) (±- - sj\ A^q) 

+ lY, [<(-g)co S [n0(k s ,q)] +A t n (-q)sm[n9(k a ,q)]] J(S; g) 
s 

+ iJ2 [A e n (q) cos[n0(k s , -q)] + <(g) sin[n0(k s , -q)]] J(S; -q)\ 



(6.19) 



where it is understood that the repeated index n is summed over. When the gauge fields are 
integrated out in the generating functional the original interaction is reproduced. Now if we 
express the currents J(S; q) in terms of the canonical boson fields we see that 



/ cos[n0(k s ,q)],/(S;g) 

q,S 

\% E<(?W^(k S) -q)]J(S;-g) 



Y E \/ Qfi s • q [a(S; q)A l n (-q) cos[n0(k s , q)] 

q,n s q>0 

+ a*(S;g)4(g)cos[nfl(k s ,-q) 



(6.20) 



and therefore 
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FL 



1 r du ^ 
V J 2^^ 



+ % «( s ;?) 

nsq>0 



Ai(-q) ( -^Ai(q) + Al(-q)(-ir { l^Ai(q) 
Qh s ■ q[A e n (-q) cos[n6(k s , q)] + A^-q) sin[n0(k s , q) 



+ a*(S-q) 



fins ■ q[A e n (q) cos[n#(k s , -q)] + <(g) sin[n#(k s , -q) 



(6.21) 



The action is now in a form which allows the canonical bose fields, a and a*, to be integrated 
out and the resulting generating functional is expressed as 



Z[tC\ = J DAiDAiexp-SdA] 



x exp - 
+ r(S;g) 



^E E 4 -^5 q - tr(S; g)C(S; -q)A-(g) + g(S; g)^(S; q)^(-g)] 



%u - v F n s ■ q 



-£(S;g) 



(6.22) 
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here the index m ranges over (£,t) and the associated 2- vectors b™(S; q) have components: 

#(S;q) = cos[n0(k s ,q)] 

£(S;q) = sinH(k s ,q)]) . (6.23) 
Using this notation the gauge action can be written as 



3g[A] 



duj 
V J 2^ 



E K (-q) KZ,{q) A™,{q) 



(6.24) 



q;n s q>0 



where the matrices K nn > are given by 



K nn'(l) = ("I) 

and 

C(?) = H) B 



Snn' __ A fj S • q 

/„ (2tt) 2 y ( iuJ - v *f^s ■ q) 



cos[n#(ks, q)] cos[n'0(k s , q)] 



(6.25) 



<W(1 - ^ 



n,0j 



A 



E 



n s ■ q 



(2tt) 2 ^ (zcj -t>£,n s ■ q; 



sin[n6 l (k s , q)] sin[n'0(k s , q)] 



(6.26) 

In these equations the q-sum is now restricted to ns ■ q > 0, and the sum over S (integral 
over angles) runs free. Recognizing that AJ2s = ^kp, equation |6.25| can be written in the 
more familiar form: 



kUq) = (-i) n 



S -fL-N*(0) 

Jn 



d6 cos 9 
2tt x — cos 6 



cost no) cos [n 



(6.27) 



and a similar equation for K^ n ,, here as usual x = u/(vp\q\) and N*(0) = m*/2n is the 
density of states of spinless fermions in two dimensions. 

The gauge fields can now be integrated out to provide an explicit expression for the 
generating functional 



z[t,e 



exp 



. r duj 
^ l 2^ 



q,nsq>0 



Er(S;?): 



1 



L S 



iuj - v F n s ■ q 



x 



Vm ^ 1 -«t ;9 ) 



(6.28) 



iuo — Dpiix • q 

Here the gauge propagator is given by the matrix inverse of the .fT-matrix: 

B m {q) = [K"^)] -1 . (6.29) 

The boson propagator now can be determined directly by differentiating the generating 
functional, equation |6J| : 

G fl (S,T;g) 



iuj — v F ris ■ q iuj — i/pns • q 



iu — v F n T ■ q 

(6.30) 
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and the fermion propagator follows immediately upon exponentiating: 



G F (S;x, r) 



A 



(27r) 2 n s ■ x + ivpT 



exp 



E 



^/.i(q.x-<JT) 

2tt 1 



bgSj -q) £ff B ,( g )6ff(S;q) 



A 



(27r) 2 n s • x + ivpT 



exp[— 5G^(S; x, r)]; xj_ = x x fi s = 



(6.31) 



This key result is exact within the bosonization scheme. The condition xj_ = can be 
replaced by |xj_|A < 1 provided we control the logarithmic infrared divergence of the free 
boson correlation function by placing the system in a large box. The weak divergence, which 
is a consequence of treating the Fermi surface as locally flat within each patch, presents 
no real difficulties in practice and is ignored in the following. For |xj_|A > 1 the fermion 
Green's function vanishes since bosons separated by large perpendicular distances, xj_, are 
uncorrelated. For convenience in the second line of equation |6.31| we have introduced the 
notation SG^S; x, r) to denote the modified Fourier transform of the additive correction to 
the free 0-boson propagator due to interactions 

«7,(S;x,t) = X^T^p / ^[e^-^ - l]^(S;q,u,) (6.32) 
where we can identify 5G«^(S; q, uj) as 



6-(S;q) D-,(g)^(S;q) 



\IUJ 



v F n s ■ q) 



(6.33) 



In equation |6.31| the momentum integral ranges over (— A/2, A/2) in each direction. By con- 
trast in the free part of the boson and fermion propagators the perpendicular momentum 
ranges over the much larger interval (—A/2, A/2). The reason for this lies in the fact that 
the Fermi surface normal points in a different direction in each patch on the Fermi sur- 
face. As interactions couple together different patches, and as the patches are squat boxes 
with dimension A D_1 x A, only wavevectors |q| < A are permitted by the geometry of the 
construction. 

To complete the program of computing the fermion two point function, the in-patch 



Green's function equation 6.31 is summed over all patches on the Fermi surface, 



Gi 



^'G f (S;x,t). 



(6.34) 



As in section 0, the prime indicates that the sum is only over patches for which |x x n s | A < 1. 



The expression for the fermion Green's function given here, equations |6.31| and |6.34|, is 



identical to that found by Castellani et al.[42| who made use of Ward identities derived in 
the fermion basis JO]] to obtain equation |6.31 . 



In general the calculation of correlation functions is complicated as it involves large ma- 
trices, but to illustrate the method, as is common practice, we can truncate the Landau 
expansion and set all parameters except f equal to zero. First, to demonstrate the nonper- 
turbative nature of the approach explicitly, we consider the propagator of a fermion liquid in 
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one spatial dimension. In one dimension there are two Fermi points and fermion quasiparti- 
cles can only propagate to the right or to the left. In this instance the generating functional 
is given by 



Z[t,?] = exp-£ f^\c(R,q) . 1 ; t(R,q)-e(L,q) , 1 ; £(L, 
g J 2n [ iu — v F q iu + v F q 

£(R,q) C(L,-q)\\ 



q) 



2ir \iu-v F q iu + v* F qj 00 W 1 
where the single element of the matrix K Q0 is given by 



iu — v* F q iu + v F q 



(6.35) 



K 00 (q) 



q 



(6.36) 



f 2n \iu — v* F q iu + v* F q / 
The boson propagator for quasiparticles propagating to the right, for example, is given by 

1 q K m \q) 



Gb(R; q) = . , , . 

iuj - v F q 2tc (iu - v F q) 2 

which after some algebra can be rewritten as 

a b 



G B (R-q) 



+ - 



where 



iuj - qv* F (l + h/KV F y/ 2 iu + gv* F (l + / /™f) 1/2 
(1 + f /2wv* F ) " 



a 
b 



1 ± 



(6.37) 



(6.38) 



(6.39) 



Now v' F = v F (l + fo/2nv F ) is the renormalized Fermi velocity due to intrapatch scattering. 
Introducing the dimensionless coupling F = f /(7tv' F ) the boson Green's function can be 
expressed as 



G B (R,q) 



cosh 2 T] 



+ 



sinh 2 7] 



where 



iu — qv' F (l — (^ i ) 2 ) 1 / 2 iu — qv' F 

1 



l-(f)2)l/2 



cosh r\ 



1 + 



(*b/2) s 



(6.40) 



(6.41) 



reflecting the fact that interpatch interactions further renormalize the speed of right and left 



moving quasiparticles, vf = v' F Jl — (F /2) 2 . The fermion propagator can now be found 
analytically using equations |672| and B]3| as 



G f (R;x,t) =ie lkpx 



IX" + VvT 



(x + ivpr) 



(6.42) 



The result is exact and identical to that found by diagonalizing the bosonic Hamiltonian 
via a Bogoliubov transformation ||. It illustrates the well-known fact that in one spatial 
dimension even short-range interactions destroy the Fermi liquid fixed point. The propagator 
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is non-analytic in the coupling constant /o, the anomalous exponent £ = sinh 2 ^, and the 



simple pole of the noninteracting Fermi gas at x 



-iv 



pT is replaced by a branch cut. The 



distribution function for fermions near the Fermi surface is obtained by Fourier transform of 
the equal-time (r = 0) propagator: 



nnik) = n^(/cp) — Csgn(A; — kp)\k — kp 



2C 



(6.43) 



where C is a nonuniversal constant. The occupancy no longer has a step discontinuity at 
the Fermi surface, rather it is continuous, although non-analytic, with a divergent slope at 
the Fermi points, characteristic of a Luttinger liquid. 

We have shown that bosonization transforms the problem of interacting fermions into a 
problem of Gaussian bosons which in one dimension can be solved exactly. In one dimension 
the conventional Fermi liquid picture breaks down and the fixed point is a Luttinger liquid. 
In the remainder of this section we use the multidimensional bosonization approach to study 
the properties of interacting fermions systems in higher dimension, especially two spatial 
dimensions. In D = 2, the zero-zero elements of the inverse gauge propagator are given by 



fo 



where 



xo(q,^j 



iV*(0) 



dO cos 9 
2w cos 9 — 



It follows from equation |6.30| that the boson propagator is given by 

^S,T 



G B (S,T;q) 



+ 



A Vns-q 



fo 



vnT-q 



iuj — v F us ■ q (27r) 2 iu — vphs 
In particular the in-patch boson propagator is: 

G B (S;g) - 1 A 



q \1 + foXo(q)J iu - v* F h T ■ q 



(6.44) 



(6.45) 



(6.46) 



n s ■ q 



fo 



iuj - f^n s ■ q 
iu - v* F h s ■ q 



(27r) 2 (iuj - v F h s 

A r , fo 
" (2ny [US ' q> 1 + f oXo (q) 



q) 2 1 + foXo(q) 



(6.47) 



Analytically continuing to real frequency, the longitudinal susceptibility is given by the 
Lindhard function: 



N*(0) 



1 - 



9(x 2 - l)\x\ M i0(l ~x 2 )\x\ 



(6.48) 



where x = u;/(f^|q|). In the so-called g-limit (the limit ui — > is taken before q — > and 
hence the ratio o;/|q| — > 0) the Lindhard function has an imaginary part: Xo ~ N*(0)[1 + 
ia;/(?/p|q|)]. The bosons have a finite lifetime because they scatter between patches. In this 
limit the boson propagator is expressed most naturally in terms of the scattering amplitude 
A Q = Fq/{1 + F ) where as before the dimensionless coupling F = N*(0)f . 



G B (S; q,uj) 



^sVi ^ - v F n s ■ q 



1 - 



A 



2nk f 



(6.49) 
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The bosons have a finite lifetime, as they scatter between patches, but unlike in one dimen- 
sion, the renormalization of the boson velocity due to interactions amongst the bosonized 
shell of excitations with energy less than VpX is insignificant, vanishing in the N = 2nkp/ A — > 
oo limit. 

To consider the effect of long range interactions, for example the Coulomb interac- 
tion V(q) = 27re 2 /|q|, we simply replace the Landau coupling f in equation |6.47| by 



the momentum-dependent interaction^] V(q). In the limit, the effective interaction 
V(q)/[1 + Xo^(q)] is screened and the analysis given above holds. In the opposite limit, the 
o;-limit, the collective excitations of the interacting fermion system are given by the zeros of 
the boson propagator where the self-energy diverges, detK = 0. In particular, xo ~ — V (2x 2 ) 
and therefore 

detK = 1 - N*(0) I ] (6.50) 




which vanishes when uj 2 = 27re 2 n/|q|/m* where rt/ = k F /4:7r is the density of spinless 
fermions. However, it is well known that the frequency of the longitudinal collective mode, 
the plasma frequency, should only depend on the bare mass. This requirement, known as 
Kohn's theorem []59|], is an immediate consequence of Galilean invariance, as large numbers 
of electrons sloshing back and forth in response to an applied oscillating electric field carry 
momentum corresponding to the product of their speed multiplied by the bare electron mass. 
To see that Kohn's theorem holds within the bosonization scheme, recall that the effective 
Hamiltonian of a fermion system interacting via long range forces describes a gas of low en- 
ergy excitations coupled by both long range Coulomb interactions and by short range Fermi 
liquid forces. For our purposes it is sufficient to set all of the Landau parameters except fi 
equal to zero. K is now a three by three matrix whose determinant is given by 



detK 



V(q) 1 + Xo, Xeo, 
Xeo, fi 1 + Xu, 

o, o, /r 1 + xu 



(6.51) 



where the susceptibilities are given by 

Xm{x) = N*(0) 



d9 cos 2 6 N*(0) 



2tt cos 9 — x 2x 

. d6 cos 3 6 3N*(0) 
rt*) = ^(0) / ^^~ x « (6.52) 



, d9cos9sm 2 9 N*(0) 
Xtt(x) = N*(0) / « -— V (6-53) 



and 

= N*(ri\ ( _ 

2tt cos 6 — x 8x 2 

the approximate expressions are given in the tu-limit. The determinant vanishes when 

^ = -Xo(x) + hxl(x) - (^) 2 (1 + y) , (6-54) 

but (1 + Fi/2) = m*/m and therefore the zeros are located at uj 2 = (2ne 2 nf /m) |q|, which is 
the correct plasma frequency of a two dimensional Fermi liquid. 
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Now we can determine the quasiparticle propagator of fermions interacting via the short- 
range Landau interaction f using: 



G F (S; x, r) = G° F (S; x, r) exp[-<5G^(S; x, r)] . 



(6.55) 



In general this is a complicated procedure; but for our purposes here it is sufficient to expand 
in powers of SG^, as the expansion converges quickly for low fermion energies. The Fourier 
transform can then be taken directly, 

G F (S; k, u) = ]T j ^ G% (S; k - q, u - v) [5^5{y) + SG^(S; q, v) + . . .] (6.56) 

fa 



where 



<5G (S;q, u) 



1 



(6.57) 



(iu - v F h s • q) 2 1 + /oXo(q) 
If we further expand in powers of /o, the first order correction is purely real and leads to 
only an infinitesimal velocity shift to the fermion propagator. The leading contribution to 
the imaginary part of the fermion self-energy appears at next order, 



SG^S; q, v) 



/o*o(q, v ) 



[iv - v F n s ■ q) 



(6.58) 



Now the only dependence of 5G$ on the perpendicular momentum q± = hs x q is through 
the susceptibility Xo(q, v): 



Xo (q,f) = JV*(0) 1 



v sgn(z/) 



(6.59) 



yj" 2 + (v* F q\\) 2 + (v* F q±Y< 
here, as usual, q» = fig ■ q. Carrying out the integral over the perpendicular momentum in 



equation 3.56, and retaining only the imaginary part we find: 



/•A/2 i\T(U) 

Im / dq ± Xo(S; q, v) = — — v sgn(z/) In 

J-X/2 V% 



N*(0) 



A 2 



(6.60) 



and therefore the contribution to the quasiparticle Green's function is 

. fi N*(0) f" f W dq\\ 



5G F (S; k\\,u 



dv — 1 zVsgn(i/) In ((v/v F ) + q\ 

J J — A/2 s-i7T 



(27r) 2 v* f Jo J-x/2 2ty 
1 1 



%{v - u) - v* F {q\\ - k\\) {iv-v* F q\\f 



(6.61) 



The integral over q» can be carried out by complex integration and the final integral over v 
evaluated exactly with the result that 



faY N*(0) sgn(u;) 



2tx J 2(v* F ) 2 (vphi - iujf 



2 ( iuj — v* F k\\ 



In 



iuj — v* F k\\ 



v* F \ 



+ 



[10J 



' iuj — v* F k\\ 



In 



iu + v F k\\ 



(2iu> - v* F k\\] 



(6.62) 
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Analytically continuing to real frequencies we find that the imaginary part of the fermion 
self-energy is given by 



ImEjr(S; k\\,u) 



7o\iv*(o) , 

2^ 2^f^^ 



J2 + (<*- V*fH 



4 



In 



u — v* F k\\ 



+ 



' uj — v* F k\\ 



. (u) + v F k\\\ uj „ , 



(6.63) 



The appearance of the combination u — v F k\\ in the self-energy is the result of linearizing the 
fermion dispersion about the Fermi energy. This form is not expected to survive in detail, 
either within the fermionic representation pDfl , or perturbatively within bosonization [ 17 , 3"8 |, 
when the quadratic term in the fermion dispersion is included. This is not surprising, as 
both the quadratic dispersion and the self-energy are equally irrelevant in the RG sense. A 
consistent calculation of the self-energy must retain the quadratic term. 

The imaginary part of the self-energy at the quasiparticle pole is the inverse of the 
quasiparticle lifetime. The location of the pole has been shifted from its bare value to 



10 



v' F k\\ due to renormalization of the Fermi velocity v' F 



4i + 



rr^-l. Note that 



FL X 1 (27rfc F ) 1+F J 

A, not A, appears in this expression because it is the patch-diagonal piece of the current- 
current interaction, (fo/V) J(S; q) J(S; — q), which is responsible for the renormalization of 
the velocity. As a result the imaginary part of the self-energy at the pole is given by 



ImSpfS; fen 



/o 2 N*(0) 



(2tt) 



2vf 



sgn(^j 



u 2 \n 



AF u 2 



\ 2 k f 



*2 



(6.64) 



the form of which we recognize immediately from previous work on two-dimensional Fermi 



liquids 0, 0, §3], §4 



This quantity is always negative since u 2 <C \ 2 k F /A. Despite the 
dependence of the self-energy on In a;, the weight of the quasiparticle pole, Z^, is unchanged 
in the |k| — » k F limit. This is as expected, since as noted above short-range Fermi liquid 
interactions make only irrelevant contributions to the quasiparticle self-energy and do not 
destroy the Fermi liquid fixed point. The result is not changed if we consider the more 
general case of nonlocal short range interactions /(S,T;q) as the additional momentum 
dependence only gives rise to irrelevant operators which do not change the behavior of the 
system. 

The question remains as to whether or not singular long-range interactions can destroy 
Fermi liquid behavior. In fact Bares and Wen|65| have shown that for a sufficiently long- 
ranged interaction, logarithmic in real space, Fermi liquid behavior in two dimensions does 
break down. Here we show how this comes about within the bosonization framework and 
then discuss why Fermi liquid theory remains intact for fermions interacting via Coulomb or 
shorter-ranged forces. As bosonization is a nonperturbative approach which does not pre- 
suppose the form of the fermion propagator it provides a particularly appropriate framework 
within which to address this question. 



|q| 



D-l 



If we consider a gas of fermions interacting via Coulomb interactions V(q) oc e 2 
then the correction to the boson Green's function is given by equation |6.57| , slightly gener- 
alized to include momentum dependence: 

V(q) 1 



5G<p(S; q,uj) 



V(q)xo(q) 



[ILO 



v* F h s ■ q) 2 



(6.65) 
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In dimension greater than one, for small q and small to, with uo <C |q| screening is complete, 
5G^ no longer depends on V(q), 

5G (S; q, to) « Xo ^ (9) . (6.66) 
{ito - vj,n s • q) 2 

The imaginary part of the fermion self-energy can be obtained in the manner we have de- 
scribed above and is given in two dimensions by 

. , 1 / \ ( 2 2ttlu 2 A lo 2 \ . . 

ImSH-) P o,e - J^MmO) 8 ^ V WM~f - T ) (6 - 67) 



and in three dimensions by 



A 



ImS F (w) po i c oc ^ 2 2 to sgn(cu) + 0(u /k 6 F ) (6.68) 



in agreement with the well known results of Fermi liquid theory |TI]. As in Fermi liquid theory 
the momentum cutoff appears in this result. In the random phase approximation (RPA) the 
cutoff A is replaced by an energy of the order of the plasma frequency p5[]. Azimuthal non- 



planar scattering processes are included in RPA, accounting for the shorter lifetime. These 
processes are neglected in our D = 3 calculation, though it should be possible to treat them 
perturbatively. 

It is clear from this discussion that the breakdown of the Fermi liquid behavior can only 
be brought about, if at all, by high-energy processes for which screening is ineffective. To 
resolve this question we determine the equal-time fermion propagator of a gas of fermions in 
two dimensions interacting via the logarithmic potential V(q) = g/q 2 introduced by Bares 



and Wen [pql . Then 

^(S; x, r = 0) = - E (e- - D ~ ^ , 2 (6-69) 

l^x -%M 2vr (iu - v* F n s ■ q) 2 

where 

K m \q) = : - , ?{f . r (6.70) 

i + (^/q 2 )xo(q,^) 

Xo{ x ^> 1) ~ N(0)/2x 2 and therefore in the cu-limit, 

= - ,. ,f , ■ (6.71) 

{iLOf — griffm 

For this super long-range potential, the collective excitations of the system, plasmons, have 



a finite energy, given by the zeros of K 00 ,lu p = \Jg^- Here we have adjusted the density 
of states and Fermi velocity such that the plasma frequency depends on the bare mass as 
it should by Kohn's theorem. The cutoff A has to be chosen such that the energy cutoff is 
much larger than the plasmon energy v* F \ ^> u p . Obviously bosonization is less accurate for 
systems with a large plasmon gap, nevertheless we expect it to be qualitatively correct. Now 
since lj > v F \q\ we may replace (iu — v F hs • q) 2 by ui 2 in the frequency integral and carry out 
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the integrations. The frequency integral gives a constant due to the plasmon pole at finite 
frequency. It remains to evaluate the two dimensional integral over the momentum q which 
because of the factor of q 2 = q±+qf\ in the denominator gives rise to a logarithmic dependence 
on x\\ = ns • x. To obtain the fermion quasiparticle propagator SG^ is exponentiated and we 
find 



A e 



ikpx 



\\x\\ | ^ , Axy 1 
Axii < 1 



(2tt) 2 

G F (S;x h r = 0) = { (6.72) 



where the exponent ( is given by: 




(6.73) 



The nonanalytic dependence of Gp on x\\ reflects the destruction of the quasiparticle pole 
by the emission and absorption of plasmons. The power law dependence with exponent ( is 
a consequence of the nonperturbative treatment of the interaction within the bosonization 
approach. The perturbative RPA approach in the fermion basis gave instead a logarithmic 
dependence on x\\. Bares and Wen|J5| then used a variational wavefunction to argue that 



the logarithm should be exponentiated into a power law with the same exponent as obtained 
above. Multidimensional bosonization obviously allows us to obtain the power law directly. 
As we have discussed, the appearance of an anomalous power law in the Green's function 
results in the elimination of the discontinuity in the quasiparticle occupancy at the Fermi 
surface, replacing it by a continuous but nonanalytic crossover. The quasiparticle occupancy 
decreases from one to zero over a momentum scale of the order of A. The derivation of 
the current algebra only requires that the occupancy change over a momentum scale of the 
same order, providing a consistent picture. It should be noted, on the other hand, that 
consistency does not necessarily mean that the result is correct. For example, it is possible 
that the super long-ranged interaction actually induces Wigner crystallization^]. Either 
way, the destruction of the Fermi liquid is complete. 

It is instructive to repeat the calculation of the equal-time Green's function for the 
realistic case of a fermion liquid interacting via Coulomb forces in both two and three 



dimensions ||31|| . In two dimensions the plasmon frequency is proportional to J |q|. In this 



case the integral over u is no longer constant, as it was for the Bares and Wen potential, 



but rather it is proportional to J |q|; as a result the momentum integral vanishes as JX/kp. 
In three dimensions the plasmon frequency is non-zero but the integral over q is again sup- 
pressed, in this instance because of the reduced phase space at small momentum, and as 
a result dG^ vanishes as X/k F . In each case the quasiparticle propagator has the standard 
Fermi liquid form. The weight at the pole undergoes no additional renormalization beyond 
the factor of Z k due to integrating out the high-energy degrees of freedom, equation |2T2| . 

In this section we examined the effects of both short and long range longitudinal inter- 
actions on fermion liquids. We found that if the superconducting instability is suppressed 
then the Fermi liquid state is the ground state of a degenerate gas of fermions interacting via 
short-range forces, or by Coulomb forces, in either two or three dimensions. Bosonization 
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provides a means to go beyond an assumed Fermi liquid form for the quasiparticle propaga- 
tor. Indeed had we considered fermions with spin we would have found that the bosonized 
Hamiltonian separates into two parts, charge and spin, leading to the possibility that, as 
in one dimension, the quasiparticle propagator would reflect spin charge separation. Spin 
charge separation would destroy the Fermi liquid, as the key element, the existence of a pole 
in the single particle Green's function with non-zero spectral weight, would be replaced by a 
branch cut. In dimension greater than one this does not occur because in the A — ► limit the 
location of the pole in the boson propagator is unchanged from its free value. Consequently 
the spin and charge velocities are equal and both degrees of freedom propagate together in 
the usual quasiparticle form. On the other hand, the Fermi liquid form is destroyed both in 
one dimension, where the ground state is a Luttinger liquid, and in the case of the superlong 
range interaction in two dimensions studied by Bares and Wen. The non-Fermi liquid fixed 
points of these two examples are accessible via bosonization. 



7 Two Particle Properties: Density Response 

In the previous section we discussed the single particle properties of an interacting fermion 
system in some detail. Here we turn to the consideration of two-particle properties, specifi- 
cally the density response. We are concerned with two limits in particular, low-wavevectors, 
|q| ps and high-wavevectors, |q| ps 2/cp. 

In the low- wave vector limit we show that the classic RPA result for the density response 
is exact within multidimensional bosonization. In fact in early work in the fermion basis this 



result was anticipated by Geldart and Taylor |68|j. These authors showed that, in the calcula- 



tion of the density response, self-energy corrections to the fermion propagator were canceled 
by vertex corrections at leading order in an expansion in the Coulomb interaction. This 
question has been addressed most recently by Kopietz, Hermisson, and Schonhammer p7| 
who were able to show that, provided the quasiparticle dispersion is linear, this result holds 
to all orders in a loop expansion. It is not immediately obvious a priori that a bosoniza- 
tion scheme with the restriction |q| < A < A can describe processes with large momentum 
transfer. In fact, as in one dimension, multidimensional bosonization can describe processes 
with either small or large momentum transfer provided they are of sufficiently low energy, 
| a; | <S VpX. To illustrate this point we determine the density response of a Fermi gas at 
|q| ps 2k F which, due to the sharp reduction in the phase space for low energy scattering, is 
nonanalytic, as was first pointed out by KohnpOf . 

To begin we write the fermion density-density correlation function in the bosonic repre- 
sentation. The fermion density fluctuation is given by 

5p(x,r)= :^(x,r)^(x,r) : = £ : ^(S; x, r)^(T; x, r) : (7.1) 

ST 

where the colons denote normal ordering which in this case is equivalent to subtracting the 
uniform background density. Then the density- density correlation function can be written 

as 

(5p(x,T)5p(O,0)) = -^£(J(S;x,r)J(T;0,0)> 

V ST 
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+ (^) 2 T,^ e^ kT > x (exp^^[0(S;x,r)-0(T;x,r)] 



x expj^[0(U;O,O)-0(V;O,O)]J) (7.2) 

For small momentum transfer |q| < A <C kp only the first term is relevant, whereas the 
second term is important at large momentum transfer. We separate the analysis of these 
two regimes in the following subsections. 



7.1 Response at Small Wavevectors 

Consider the density response of a free Fermi gas at long wavelength and low frequency. For 
free fermions, the currents in different patches are uncorrelated and the density response is 
given by 

n°(q,u;) = -£(J(S;q,u;)J(S;-q,-u;)> (7.3) 



where J(S; q, u>) is given by equation [2.20| and therefore 



E(J(S;q.^) J(S; -q, -w)> = fi£( n s-q) ( a (S;q)o*(S;q)) 
s s 



E 



n s ■ q 



g iuj - v F n s ■ q 

If we consider the specific case of two dimensions and recall that v , ^ 2 ^~ — v</m r •><> 



then 



d6 cos 9 



n°(q,cu) = N*(0) / 
where again x = u/(v F \q\) and therefore 



(7.4) 
(7.5) 



n (q,u;)=iV*(0) 



\x\ 



Vx^l 



(7.6) 



which is the well known result for the density response of the free electron gas in the limit 
of small momentum. 

Generalizing to include a long range Coulomb interaction, currents in different patches 
are now correlated and the correlation function written out in full reads 



(J(S;g)J(T;-g)} = ^(n s • q)(n T ■ q) x 

+ (a(S;q)a\T;q))e(n s - q )e(n T - q ) 

+ (at(S; -g)a(T; -q))6(-h s ■ q)0(-n T ■ q 

+ (at(S; -g)ot(T; g))0(-n s • q)fl(n T • q 



(a(S;g)a(T;-g))e(n s -q)0(-n T -q) 



(7.7) 
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and contains contributions from all possible two point boson Green's functions. All of these 
functions can be obtained from the generating functional, which in this instance is given by 



n s ■ q 



A(-q)A(q) 



x 



expj] 



n s q>0 



iQy/hs ■ q 



- T(S;g)- 



zcj - ^n s • q 
£(S;g) 



[A(-g)£(S;g)+A(g)f(S;g)] 



(7.8) 



zo; - t£n s • q 

here A(g) is the longitudinal gauge field introduced to decouple the Coulomb interaction. 
The correlation functions in equation \l.7\ are found in the usual way, for example 

S 2 In Z 



(a(S;q)a\T;q)) 



5e(S;g)^(S;g) 



fi,/(n s ■ q)(n T ■ q) 



(iuj - v* F h s ■ q)(iuj - v* F h T ■ q^ 
Likewise, the anomalous propagator is given by 



(A(q)A(-q)) 



(7.9) 



(a(S;g)a(T;-g)) 



£lj(h s • q)(n T • q) 



(iu - v* F h s ■ q)(ia> - v* F ii T ■ q 



(A(q)A(-q)). 



(7.10) 



Combining all of the terms in equation |7J we obtain 

fin s • q 



(J(S;q,u) J(T;-q,-u;)) 



+ tt 



sD-l 

^ - ^n s ■ q ' 

A (n s -q)(n T -q) 



x 



(27r) 2 {iuj - v F n s ■ q)(iw - v* F n T ■ q) 
V(q) 



Then summing over S and T we find 

n(q,w) 



1 - V(q)II°(q) 



n°(q,q;) 
l-y(q)n°(q,u;) 



(7.11) 



(7.12) 



which is identical to the random phase approximation result for the density response, and is 
exact in the limit q, u — > 0. 



7.2 Response at |q| « 

We have seen that bosonization successfully captures the physics of the long wavelength 
response of an interacting fermion system. As a further test of the formalism we consider 
the density response at short wavelengths, in particular around |q| = 2k F . At wavevectors 
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Figure 10: Wavevector of magnitude 2kp just spans the Fermi surface. 



greater than the spanning vector of the Fermi surface it is no longer possible to excite particle 
hole pairs at zero frequency; see figure |T0[ Consequently the density response is nonanalytic 
at |q| = 2k F , as was first pointed out by Kohn|69[ . 

At first sight it is not obvious how large momentum transfer processes can be described 
within a bosonization framework with the restriction |q| < A < kp. The reason why this is 
possible is made clear on noticing that scattering between patches is allowed, and that in 
fact the second term in equation |7]2| contains these processes explicitly; recall that anoma- 
lies in 2kp scattering are described correctly by standard one dimensional bosonization [^SJ. 
Bosonization in one and higher dimensions can describe processes with either small or large 
momentum provided they are of low enough energy to lie inside the shell depicted in figure 



Since for free fermions the bosons in different patches are independent, the correlation 



function in the second term of equation [7.2| is zero except when V = S and U = T and the 
contribution to the density response function is given by 



n°fx,r) 



XVaJ 



\ e i(ks-k T )-x^ e i^[ v iS:x.r.- 



A 



D-l 



(2< 



D 



E 



n s • x + ivpr 



A 



D-l 



(2vr; 



D 



E 



(T;x,r)] e i^f[<«T;0,0)-</,(SAO)] 



n T • x + iv F r) 



(7.13) 



Implicit in the last line is the restriction |xj_|A < 1. As we are concerned with the response 
at |q| pa 2kp, only patches which are nearly opposite to each other on either side of the Fermi 
surface participate in the scattering, fig pa — n T . Therefore we may carry out the sum over 
patches in each of the two terms in equation |7.13| keeping in mind that the spatial directions 
perpendicular to the Fermi surface normal are the same for both terms. To simplify the 
calculation we consider the static response, u — 0, only and specialize to the case of a 
circular Fermi surface in two dimensions. Integrating the density response over time r we 
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obtain 



dd X f d6 2 e ikFr(co B 6 1 -co S e 2 ) 



27r r(cos#i — cos 82) 



[sgn(cos6'i) — sgn(cos^2)] (7.14) 



where 8\^ are the angles between x and ks, and x and kr, respectively, and r = |x|. Now 
we can focus on the case of |q| ~ 2kp by limiting the range of integration over the angular 
variables 8^2 to angles near 8\ = and 82 = vr. In this limited range the cosines appearing in 
equation |7.14| may be expanded in a Taylor series as cos 6*1,2 = ±(1 — 8\ 2 /2) and the Gaussian 



integral over angles can be carried out. Finally take the Fourier transform to q-space. Then 
if 8 is the angle between q and x, the density response function is given by 



n°( q ,^ = o) 



[N*(0)f 
2m* 



rdr 



2lT , / .expi(2/c i r — Iql cos( 
do 



Now 



- r d8e~ l ^ rcose = Jo(|q|r) 

7T JO 



— 1 — — cosl o r ) 

7r|q|r 4 



|q|r ^> 1 



and therefore 



n 2 fc F (q,^ 



0) 




00 exp[(i(2k F - |q|)r] 
dr Z3/2 



(7.15) 



(7.16) 



(7.17) 



Here we have assumed |q| ~ 2kF and we have also introduced the ultraviolet cutoff a to 
eliminate the artificial divergence at small r. Finally, asymptotic analysis of the radial 
integral gives the nonanalytic contribution to the density response function 



e i(2fc F -|q|)r 



.3/2 



dr « — 2\fniJ |q| — 2k 



and therefore 



n 2 fc F (q,w = 0) 



V2n 



■Itt 



|q|-2fc F , > Y 
2k F h °> 



C. 



> 2k F 
< 2k F 



(7.18) 



(7.19) 



Here C is a constant which can be determined by a more careful evaluation of the integrals 
in the small r region. Continuity across |q| = 2kp ensures that the same constant appears 
in both limits of equation |7.19| . The Kohn anomaly of equation |7.19| agrees precisely with 
that found directly in the fermion basis |7(J. Repeating the calculation in D = 3 we show 
easily that the static response is given by 



n 2 fc F (q,cu = 0) 



k F m* 2k F 




|q| 


8tt 2 


q| 





In 



|q| — 2k F 



2k, 



+ 



(7.20) 



Which again agrees with the exact result^] found in the fermion basis in the limit |q| — > 2kp. 
The Kohn anomaly persists even if longitudinal interactions are included. This has already 



been anticipated in equations |7.19| and |7.20| where the effective mass appears rather than 
the bare mass. 
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8 The Half-Filled Landau Level 



The two-dimensional electron gas in a magnetic field with even denominator Landau level 
filling fraction has been the subject of extensive study. For an overview see [71]]. Experiments 



at even denominator filling point to a compressible metallic state |72|, [73], |74|, [75|, [TJ| rather 
than the incompressible state found at odd denominator fractional filling. Furthermore, 
the remarkable similarity between phenomena observed at high magnetic fields near filling 
fraction v = 1/2, and those found for electrons in weak magnetic fields can be explained 
within the composite fermion framework |77||. A composite fermion is formed by attaching 



an even number of "statistical" flux quanta to an electron. Specifically, in the case of the 
half-filled Landau level, two such flux quanta are attached. Then if the statistical magnetic 
field due to each attached flux line is spatially averaged, the external magnetic field is 
effectively canceled and the composite fermions travel in zero net field. The mean-field 
ground state is, therefore, a filled Fermi sea with a Fermi momentum kp larger than that of 
the free, unpolarized, Fermi electron gas by a factor J78J of y/2 consistent with experimental 
observation ||76||. 

This picture has been formalized as a Chern-Simons gauge theory by Lopez and Fradkin[79] , 



by Kalmeyer and Zhang [pO[| , and by Halperin, Lee, and Read|78]] (HLR). In this approach 
the electrons are subjected to a singular gauge transformation converting them into com- 
posite fermions which interact with each other via a statistical transverse gauge field, the 
Chern-Simons field, in addition to the usual Coulomb or short-range interaction, and the 
mean-field theory is constructed. However the effect of dynamical fluctuations of the gauge 
field (equivalent to fluctuations in the electron density) must be addressed. In this sec- 
tion and in the following section we employ multidimensional bosonization in an attempt to 
answer this question. 

The Hamiltonian for a two dimensional electron gas of density nj in a perpendicular 
magnetic field B is given by 

H[4, c e] = ^J d 2 x 4(x) + eA(x)J c e (x) + H Cox]lomh [cl, c e ] . (8.1) 

We have dropped the spin index as we assume that the spin of all of the electrons has 
been aligned parallel to the external magnetic field by the Zeeman interaction. The vector 
potential A(x) is due to the uniform external magnetic field B which points in the z direction 
normal to the plane and m is the band mass of the electrons. We use units where c = 1 and 
the electron charge is — e. At a value B = A-KUf/e, half of the states in the lowest Landau 
level are filled. The composite-fermion theory of the half filled Landau level makes use of a 
local gauge transformation to describe the system as a collection of quasiparticles that obey 
Fermi statistics: 



c^(x) = 4(x) exp 



-i<p J rf 2 yarg(x-y)4(y)c e (y) 



(8.2) 



When is an even integer the quasiparticles are fermions. Each quasiparticle is a composite 
object comprising the physical electron together with a flux tube of integer quanta. The 
transformed Hamiltonian is given by 

H [c\ c] = f d 2 x c f (x) [- + eA(x) - a(x)] 2 c(x) (8.3) 
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where 

a(x) = 4> fd 2 y ^ (X ~ y) ct(y)c(y) (8.4) 
J |x-y| 2 

and 

V x a(x) = 2tt n/(x) (8.5) 

The average field strength of the flux tubes can now be adjusted to cancel out the external 
magnetic field. At half filling, setting = 2 



V x [eA(x) - a(x)] =e[B ^ = 0. (8.6) 




Consequently at the mean field level, the quasiparticles behave as if they are free fermions 
in zero net field; the infinite Landau degeneracy is lifted and presumably the state is stable 
against two body Coulomb interactions. The important question to ask now is whether 



gauge fluctuations or, equivalently quasiparticle density fluctuations (equation |S.5| ) modify 
or destroy the mean field state. In the following we assume that despite strong gauge fluc- 
tuations, multidimensional bosonization is applicable. Later we verify that this assumption 
is correct for the case of Coulomb or longer-range interactions. 

The bare action for composite fermions interacting via a Coulomb interaction and subject 
to the constraint, equation |B.5| , can be written 

d 



S[c*, c, Ap, a] = J d 2 x dr |c*(x) (— iA (x)j c(x) 



V 

h eA(x) — a(x) 



27T0 



+ ^/ArfVr[V x xa(x 1 rMx-y)[V y xa(y,r)] (8.7) 

As usual we work in the Coulomb gauge V • a = 0. The next to the last term in equation 
is the Chern-Simons term introduced by Lopez and Fradkin[[fj|], by Kalmeyer and Zhang pO 



and by Halperin, Lee, and Read|]78[ (HLR) so that the integration over A enforces the con- 



straint, equation |8.5| . The last term in the action is the longitudinal Coulomb interaction 
which has been written in terms of the gauge field by making use of the constraint. It is 
important to notice that in the Chern-Simons gauge theory the longitudinal interaction plays 
an important role in suppressing fermion density fluctuations. Since density fluctuations are 
equivalent to gauge fluctuations in a Chern-Simons theory, repulsive long-range Coulomb 
interactions moderate and control the gauge fluctuations. As shown below, gauge fluctua- 
tions are stronger for shorter-range interactions, and in fact completely obliterate the Fermi 
surface. In this case multidimensional bosonization is not useful. 

As before we integrate out the high-energy fermion degrees of freedom using the renor- 
malization group. We also make the replacement a — eA — > A, reflecting the fact that now 
we are studying fluctuations about zero net mean magnetic field. The resultant effective 
action reads: 



S eff [ijj*, if), A^] = d x dr <i/j*(x) [- iA (x) ) i>(x) 
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+ ^*(x) (y - A(x))%(x) + (V x A(x))| 

+ \j^- 2 J d ' x d "y dr t V * x A ( x ' r)\ 7(x - y) [V y x A(y, r)] (.8.8) 

At filling fraction i/ = 1/2 we set = 2 to cancel out the external field. As stated above, 
A should be interpreted as the fluctuation about the mean field value. The effective action 
can now be expressed in terms of Abelian charge currents 



J(S;q,u) 



as 



where 



Q\h s ■ q| [a(S; q, Ld)6(n s ■ q) + a*(S; -q, u)6(-n s ■ q)] 



(8.9) 
(8.10) 



q n s q>0 



A(g)-A(-g)} , 



^o(-g) + -4 • A(-g) 



2m* 



(8.11) 



and the bare gauge action can be read off from equation pTB 

1 



SUA,} 



1 



T/ E / r> 



V 



27T { 47T0 



[A (-?)[qxA( g )]-A„( 9 )[qxA(-g)]] 



(8.12) 



Fermi liquid interactions should also be included in the action^]: If all Landau coefficients 
except fi are set equal to zero they are described by 



Sfl = E / f£ AS; g) (ks • k T ) J(T; -g) 



do; 



(8.13) 



"F S,T;q- 

which is made gauge covariant by the replacement 



(8.14) 



As stated, for computational convenience we work in the Coulomb gauge which eliminates 
the longitudinal component of the vector gauge field, only the transverse component Ap(g) 
q x A(g) remains. 

Now we can determine the composite quasiparticle two point Green's function. It is 
important to note that the propagator for physical electrons is not the same as the propagator 
for the composite particles, as the physical electrons do not have attached flux [[78| . The first 
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step as usual is to determine the boson propagator; for the moment we set the short-ranged 
part of the interaction, the Landau parameters generated upon integrating out the high- 
energy fields, to zero. We construct the generating functional and integrate out the bose 
fields 



Z[tC\ = jDA,ex V -S G [A^ 



x exp 



Q n s q>0 



V[icu - v* F h s ■ q 
+C(S;q)[iA (q)+v F sm{9)A T (q)} 

ex P -/^ E ^(S;g) 



[Z(S;q)[iA (-q)-v* F sm(9)A T (-q)) 



q,n s q>0 



iuo - v F n s ■ q 



<(S;q) 



(8.15) 



where the gauge action 
S G [A,] = S G [A, 



2^ 



E E 

q n s -q>0 



[iA (-q)-v F sm(e)A T (-q)}- 



tth s ■ q 



iuj - v pn s ■ q 



x [iA o (q) + v F sm(0)A T (q)] 



n t 



2m* 



A T {-q)A T {q) 



(8.16) 



In writing equations |8.15| and |3.16| we have used the fact that 



k s ■ A(q) 



rrr 



-^qx A(q) sinfl(k s ,q) 

m* 

v* F sin(0) A T (q) 



117) 



where as before #(ks, q) is the angle between ks and q. Now in equation |8.16| restrict the q 
integral and let the S sum run free; then we can write the gauge action as 



Sg[A»] = Ij: £ J d £A,{-q)K^{ q )AM 

S nsq>0 

where the matrix K is given by 

K(,)=l * »■<«)- 



(8.18) 




(8.19) 



and \x and v are summed over (0, T). Here xo is the longitudinal susceptibility given by 
equation |6.48| and xt, the transverse susceptibility, is given by 

d6 sin 2 6 cos rip 



xnq,wj 



v*/N*(0) 



2tt cos 6* — x 



m 



-^ 2 iV*(0)[x 2 - \x\Vx*^19(x 2 - 1) - i\x\VT^x?9(l - a; 2 )] ; (8.20) 



as before, x = u/(v F \q\). Integrating out the gauge fields and differentiating the logarithm 
of the generating functional with respect to the source fields £ and £* we find the boson 
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propagator 

Gb[s, T; q) — — 4t— + E >:(S; q)fl""( ? )^(T; q) 



iuj n - u^n s q p^tu - v* F n s ■ q zo; — v* F n T ■ q 

(8-21) 

where the vector 6^(S; q) has components 

6 (S;q) = 1 

6 T (S;q) = z^sin0(k s ,q) (8.22) 

and the gauge propagator is given by 

D Vfi = (A u {-q)A^{q)) = (K~ l y». (8.23) 

In the g-limit, u <C fj?|q|, the gauge propagator is dominated by its transverse (T) 
component and 

&;(S; q)^( 9 )6,(S; q) « A , . (8.24) 



(27T</>) 2 / (27T(/)) 2 

In this limit the gauge propagator has an imaginary part due to quasiparticle damping. For 
the Coulomb interaction, after the analytic continuation iu — > u, 

fc*(S; q^q^S; q) » vf % , ■■ 1 (8.25) 

q ^y-xlql 

where following the notation of HLR we introduce 7 = n f / itv F N (ft) , and x = 2ire 2 / (47r 2 e<j) 2 ), 
where e is the dielectric constant. The pole is located at u = i(e 2 /£/ci?)|q| 2 . If, on the other 
hand, V(q) is a short range interaction the pole is at u oc iv F \q\ 3 k F . The location of the 
pole in both instances agrees with the RPA result of HLR[[7£|. In the opposite limit cu-limit, 
u > v* F \q\ 

u? (n f /m*)(2ir^ 2 
q 2 u 2 - {n f /m*) 2 2Ti(j) 2 

In this limit , the zeros of the denominator give the characteristic frequency of the collective 
modes of the system, the cyclotron modes, with frequency to = In^nf /m* , the cyclotron 
frequency of a free particle of mass m* in an external magnetic field of magnitude B = 
27rn/0/e. The correct result is obtained, as in Section || when the Fermi liquid interaction, 
equation |3.13 , is included in the theory. Then the zeros shift to 

\u\ = u c = 2imf(f) ^ — - + = 2nnf(j)/m (8.27) 

as predicted |59|]. The details of the proof of this result are given in Appendix |D]. 

The form of equation |8.26| is similar to the corresponding result found for the super long- 
range interaction of Bares and Wen|65|. However since a projection onto the lowest Landau 



6J(S; q)^(g)6 M (S; q) « ^ S^U^l , . (8 . 26) 
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level is implicit in the HLR theory, the cutoff must be such that v* F X <C uj c to disallow 
inter-Landau level transitions. This restriction is an auxiliary condition in the definition the 
composite fermion Green's function. Non-Fermi liquid behavior, if there is such behavior, 
arises from the form of the gauge propagator in the opposite g-limit. Although there is no 
obvious small parameter we can estimate the fermion quasiparticle self-energy from equation 
B.24| by expanding the fermion Green's function in powers of D^ v . Given this caveat, in the 
case of the Coulomb interaction, at first order in the expansion, the self-energy has the 



marginal Fermi liquid form in agreement with the RPA calculation of HLR[78 



As the expansion in powers of the boson propagator is unreliable instead we calculate 
the fermion Green's function directly from the key bosonization formula, equation [2.21| . The 
fermion Green's function equation |6.31 is a product of the free propagator and a multiplica- 
tive factor which is the exponential of the correction to the free 0-boson propagator 



G F (S; x, r) = G° F (S; x, r) exp[-5G^(S; x, r)] 



where 



5GJS: x, t) 



E 



duj 

a?' 



3 i(q-x-wr) 



1) 



6:(S;q)D^(g)6 M (S;q) 
{iuj - Vp&s ■ q) 2 



(8.28) 



1.29) 



Setting x± = 0, the integral over q\\ can be evaluated by complex integration. In the g-limit 
we find that 



5GJS; x, t) 



Fill 



2nx Jo 2tt 



00 d,U! -^-sgn(a;||)(x||+TO* r) 

e V 



In oo 



In uj 



10 



5.30) 



The final integration over frequency can be carried out using the integration formula 

/ e- f * B hi(p + x)dx=-\hiP-e> tf} E i (-iiP)], (8.31) 
Jo n 

where Ei(x) is the exponential integral function. The result, after continuation to real time, 
it - 



Xii 



-t, is 
<5G (S; x,t) 
Here 

is a momentum scale, 



vU 



ln|(x|,-^)5|-e (:c ll-^' )<5sgn(x ii ) E i [-(x,|-^)(5sgn(x|,)]+7 



5 = 



A 2 



Aev* F 4> 2 k F 



ev, 



2vre 2 



(8.32) 



(8.33) 



(8.34) 



is a dimensionless number and 7 ~ 0.577216 is Euler's constant. Making use of the asymp- 
totic expansion of the exponential integral function SG^S; x, £) reduces to 



5G (S;x,i) 



C|x||<5|ln(|x||-^t|(5), 
-C-^-M\x l{ -v F t\5), 



\X\\ 



•£'11 



v* F t\S < 1 
v* F t\6 > 1 



(8.35) 
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For simplicity, first we consider the equal-time propagator which at long distances |fig • 
x| 3> is given by 

G F (S; x) = |^ (|ns ■ x\S)~< (8.36) 

and at short distances |ns ■ x| <C <5 _1 by the usual Fermi liquid form 

G ^ = ^r^ (8 ' 37) 

At momentum scales less than 5 the quasiparticle pole is destroyed by transverse gauge 
interactions, suggesting that the v = 1/2 system is controlled by a non- Fermi liquid fixed 
point. Despite the appearance of the anomalous exponent (, the composite quasiparticle 
Green's function differs in an important way from the Green's function found in Section |6] 
for the super long-range interaction. It can be inferred from equation |8.36| that the composite 



quasiparticle occupancy drops abruptly from one to zero over a momentum scale 5 oc \ 2 /k F 
which is much less than the cut off A, as A — > 0. The basic consistency condition for the 
bosonization scheme is well satisfied in this instance. 

The carrier density of a typical Ga As / Al x Gai_ a . As sample is of the order of rif pa 
10 n cm -2 , the electron band mass mj pa 0.068m e , and the dielectric constant is given by 
e pa 12.8. The effective mass of the composite fermion quasiparticle has to be estimated. 
Since the electron states in the lowest Landau level are degenerate, the electron band mass 
does not enter. Instead the energy scale is determined by the Coulomb energy at the magnetic 



length scale|78|, |8"2 |. From dimensional analysis the effective mass of the quasiparticle must 
be therefore 

m* = (8-38) 

where C is a dimensionless constant to be determined either theoretically or from experiment. 
Experiments at v — 1/2 suggest [73. 74] that m* pa lOm^ and hence C ~ 10. This value agrees 



well with the theoretical estimate of Shankar and Murthy[^]. For values in this range we 
find C = C(p 2 /2tt ~ 0.05, quite small and possibly difficult to measure. For the purpose of 
estimating 5, we may take the cutoff A to be of the order of the Fermi momentum, A ~ kp, 
and find 5 = A 2 /{AC(j) 2 kp) ~ 0{kp) as expected since kp is the only physical momentum 
scale in the problem. Notice that the Fermi energy ep — kp/(2m*) is still much smaller than 
the cyclotron energy u c = k 2 F 4>/ (2m&) and hence the requirement that v* F X u c is respected 
even for A ~ kp. 

Let us consider briefly longitudinal interactions other than Coulomb. Because longitu- 
dinal interactions control fermion density fluctuations and hence, by the constraint, gauge 
fluctuations, we expect repulsive interactions that have a longer range than Coulomb to sta- 
bilize the Fermi liquid fixed point. Drawing from the RG approach of Nayak and Wilczek|^3| 
and Chakravarty, Norton, and Syljuasen|84j] we set V(q) = glql 4 ' -1 , as y = (the Coulomb 



case) is marginal. Repeating the calculation of the equal-time Green's function for y ^ 0, 
we find 

8G,(S; x) = [ ^ r F ^ ( X -l\ eWI - . (8.39) 

0V ' ' J 2tt Jo 2vr U* ) - y \q ± 2 +y K ' 
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When y ^0, the correction to the boson propagator scales as[^, ^ 

^(S;x)~(|x|||50 2//(2+s/) , (8.40) 

where 

6' = (x/v* Fl )(v* F / X ) {2+y)/y ■ (8-41) 

The same result has been obtained in the Eikonal approach of Khveshchenko and Stamp [jS6| . 
For y < 0, the gauge interaction is irrelevant, SG^ — > as — > oo and the Landau fixed 
point is recovered. For y > 0, however, there are pronounced deviations from the Landau 
fixed point when {SG^l ^> 1 or, equivalently, when 3> The small momentum 

regime is controlled by a non-Fermi-liquid fixed point about which bosonization can tell us 
little as the deviations from Fermi-liquid behavior are too large, much larger than in the case 
of the the longitudinal interaction of Bares and Wen discussed in the previous Section. For 
example, the discontinuity in the quasiparticle occupancy at the Fermi surface is completely 
eliminated and the formalism is not internally self consistent. The problem of a pure gauge 
interaction in two dimensions, considered by Lee and Nagaosa [jB7H and Blok and Monien|p8| 
as a model for the cuprate superconductors, is the special case y = 1. Again, this system is 
outside the scope of multidimensional bosonization. 

An important feature of the bosonization procedure is the linearization of the free fermion 
spectrum in the vicinity of the Fermi surface: ek s + q — £k s ~ Vg • q. Nonlinear terms in the dis- 
persion due to curvature of the Fermi surface are neglected and within bosonization can only 
be accounted for perturbatively. This approximation may be questioned, particularly in the 
case of a mediating transverse gauge interaction, when the interaction vertex draws its most 
important contribution from q±, the momentum perpendicular to ks. To see this, consider 
the transverse gauge propagator D^(q) ~ ry|u;|/|q| — the perpendicular momentum 

9 2 

scales as qj_ oc |u;| 2 +y and therefore it, rather than the pole familiar from the linearized ap- 
proximation, q\\ = to /v F , dominates the physics of quasiparticle propagation for y > 0. The 
question of the reliability of multidimensional bosonization has been addressed|$3| in the 
fermion basis by the use of the Ward identity approach derived by Castellani, Di Castro and 
Metzner[i2|. As mentioned in Section [2] the Green's function derived using the Ward identi- 



ties is identical to that found by bosonization provided that the relationship between density 
and current vertices, A(S) = VgA°(S), which becomes exact in one dimension, remains ac- 
curate in dimensions greater than one. In this context, we can ask whether deviations from 
this relationship are so strong that they destroy the form of the propagator found within the 
multidimensional bosonization approach. For the Chern-Simons gauge theory with y < 0, 
multidimensional bosonization is internally self-consistent as the quadratic parts of the dis- 
persion cancel out and the vertex corrections are subleading||33||. The marginal case of the 
Coulomb interaction, y — 0, is delicate and must be treated with special care. Whether or 
not multidimensional bosonization is accurate for y = is a question we address below. 

We return now to an examination of the single particle properties of composite fermions, 
in particular the physically interesting case of the marginal Coulomb interaction, y = 0. The 
starting point is the single particle propagator in real time which is obtained from equa- 



tion |8.35| . In particular we show that the Green's function obtained from multidimensional 



bosonization is similar, but not identical, to the marginal Fermi liquid (MFL) form|^|. For 



57 



simplicity we set v* F = 1 and denote x\\ by x; we also drop the prefactor of Ae ?ks ' a: /(27r) 2 . 
Then the Green's function in the two limits \x — t\8 ^> 1 and \x — t\8 <C 1 can be written as 



G F (S;x,t) 



x—t+ia sgn(x) 



exp {£|x|51n |x — t|5} ; |x — t|5 <C 1 



x—t+ia 



i exD r_f sin|g-tin . 

i sgn(x) " \ ' a:— t J ' 



(8.42) 



|s - £|5 > 1 



Ideally we would Fourier transform the Green's function directly into momentum and fre- 
quency space to extract the self-energy. Unfortunately, this task is quite difficult, so we begin 
with the observation that the pole of the free fermion Green's function at x — t has been elim- 
inated by the transverse interactions, because as x — > t, GV(S; x, t) — > \x — t\^ x ^/(x — t) — >• 0. 
To find where the spectral weight is concentrated, we expand the Green's function in the 
limit of long time: \t\ ^> \x + £xln(|x — t\8). Since in the long-time limit \x/ (x — 1)\ <C 1 we 
may expand the exponential appearing in the second line of equation |8.42| in Taylor series; 
keeping the leading terms we obtain 



GV(S;x,t) 



1 



1 



x-t + iasgn(x) 1 + ^ ln(|x - t\)6) 
1 

(x — t) + (x \n(\x — t\8) + ia sgn(x) 



(8.43) 



where the ellipses denote higher order terms that are negligible in the \t\ 3> x limit. We 
may now Fourier transform the leading term to obtain the self-energy. If we make the 
approximation that near the pole in equation |S.43| , ln(\x — t\)8) ~ ln(|x|<5) the integral over 
time can be done and 



Gp(S; k, u) 



dxdt- 



i{u)t— k-x) 



x — t) + (x ln(|x — t\8) + ia sgn(x) 
2m J dx[6{-x)6{uj) - 6{x)6{-uj)]e l ^- k)x exp[t(uln\x\8] 

\i(ux ]n{\x\8)] n 

n=0 



/oo 
dx[6{-x)0(u) - 6(x)6(-uj)]e l ^- k)x ]T 
~ — n 



rv. 



. (8.44) 



Now each term coming from the Taylor-series expansion of the exponential can be integrated 
separately and the series resummed. It is readily shown that 



dx e^-^xlndxl 



(ik — rj) 



T {n\ ln"[fo - ik)/8] + Oiln^Kr) - l k)\5})} (8.45) 



where rj > is a small term which ensures convergence. Up to subleading corrections, the 
Green's function is 



G F (S;k,u) 



2tt ~ 



u) — k 
2tt 



ln[i(u; — k)sgn(uj) / 8] r 



(a? — k) — (uj\n[i(uj — k)sgn(u) / 8] 



(8.46) 
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and hence the leading contribution to the quasiparticle self-energy is 



ITT 

E(S; k, u) = lnflw - k\/6) + sgn[(a; - k)u] . (8.47) 

At the pole k = u — (culn(\uj\/8),w k and to logarithmic accuracy the self-energy is given 
by 

£(S;k,cj) = Cu\n{\u\/v F 8) - i^Qu (8.48) 

the marginal Fermi liquid formf89|. In equation |S.48| we have restored v* F . It should be 
emphasized that v* F is the Fermi velocity obtained by integrating out the high energy degrees 
of freedom; but in a marginal Fermi liquid such as that found here there is no definite velocity; 
rather the velocity is a logarithmic function of the energy scale and is of the order of, but 
not equal to, v* F . At half filling = 2 and £ = 2v F /ne 2 . Upon setting A = k F for the 
purpose of estimating the size of the self-energy we find v* F 5 = 0(e F ), the Fermi energy. It 
should be emphasized that this result is incomplete, as the equal-time composite fermion 



Green's function, equation |8.36| , exhibits an anomalous exponent. This equal-time behavior 
is subleading, however, as the Green's function decays more rapidly as a function of increasing 
spatial separation than at large temporal separation. The approximate self-energy calculated 



above, equation |8.48| , captures only the long-time behavior of the Green's function. 



Equation |3.48| is identical to the self-energy found by Altshuler, Ioffe, and Millis via a 



perturbative loop expansion pOfl. These authors argue that the multi-loop contributions to 



the self-energy only contribute sub-leading corrections to the one-loop result. This claim 
is clearly at variance with the anomalous power law decay of the equal-time propagator 
found within multidimensional bosonization, equation |3.36| . However, a recent calculation 
by Kopietz|JT| shows that the two- loop contribution to the self-energy, S 2 (e), is actually more 



singular than the one- loop self-energy £i(e) ~ e 2 /( 2 +y) ; casting doubt on the reliability of 
perturbative loop expansions. Specifically, Kopietz evaluates the two-loop Feynman diagram 
(see figure [11]), employing bare fermion propagators with the full non-linear quasiparticle 
dispersion: 

V (p - q; -q) G (P - q, e - u) D TT (q, u) V x (p, e; q, u) . (8.49) 

here Go is the bare fermion Green's function given, setting v F = k F = / y = x = ^- f° r 
simplicity, by 

G (P, V) = , ! * pr • (8.50) 

iv — 2 (p — 1) + if] sgn^z/J 
Dtt is the transverse gauge propagator which does not renormalize significantly^, ^] and 



is therefore given by the formula obtained above using multidimensional bosonization: 

DttM= 1 — . (8.51) 

|q| 1^1 

V\ is the one- loop correction to the bare current vertex V . Unlike Vi, V does not depend 
on the external frequencies, and is given simply by: 

%;q) = qxp. (8.52) 
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Kopietz finds that S 2 (e) ~ e 2 /( 2 +y) ln(e) which is bigger than E^e) by the logarithmic factor 
ln(e). This result contradicts equation 8 of Reference |)0 . 



To check the result, the one- loop correction to the current vertex has been computed 
numerically f94}| on lattices up to size 1024 x 1024 x 1024: 



Vi(p,e; q,uj) 



duj'd 2 q' 



Go(p - q - q',e 



Vo(p-q;qO 



DM4, V (p - q - q'; -q') G (p - q', e 



10 



5.53) 



An important test of the numerical calculation can be performed by invoking the Ward 
identity. For external momentum q = 0, the Ward identity (which reflects the conservation 
of the total number of particles J95, fl"2"l ) relates the one- loop vertex correction V% to the 
quasiparticle Green's function G\ dressed with the one-loop self-energy £i(e) ~ e 2 /( 2 +f); 



Vi(-p,e;0,u) 



P,ej 



UJ 



154) 



The numerical results agree with this exact relation for < y < 1, demonstrating their reli- 
ability. The vertex can also be evaluated in the kinematic region relevant for the calculation 
of the two-loop contribution of the self-energy when the fermion propagators are dressed at 
one-loop order, £ 



(i). 



P,ej 



dujcPq 
(2tt) 3 



Vb(p-q;-q) G^p-^e-u) D TT (q,uj) Vi(p,e;q,u;). (8.55) 



For the case of very strong gauge fluctuations, y — 1, the important region is given by 
q\\ ~ lo 2 ^ and q± ~ uo 1 ^, the vertex correction is finite and nearly constant as a function of 
the gauge frequency u: 



Vi(p,e; 9|| ~^ 2/3 ,g ± ~ uj 1/3 ,uj) 



const. 
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If this were the whole story, it would be a self-consistent demonstration that the leading 
behavior of the full self-energy is given simply by the one-loop form, as stated by Altshuler, 
Ioffe and Millis|9(J and by Stern and Halperin[|96|. The Fermi liquid fixed point would be 
destroyed, the quasiparticle weight would vanish, and a type of Migdal theorem^, [13] would 



hold as vertex corrections would make only subleading contributions. However, unlike the 
usual Migdal theorem where the electron mass divided by the ion mass is the small parameter 
which limits the size of the vertex corrections, in the gauge theory there is no obvious small 
parameter. 



Kopietz [91] points out, however, that in the kinematic region relevant for bare fermion 
propagators, qu ~ u and q± ~ cu 1 ^ 2 , there is a contribution to the one-loop vertex correction 
V\ which, in contrast to the bare vertex Vq, changes sign under: 

e -> -e (8.57) 



or 



UJ — > —UJ . 



(8.58) 
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p-q 




V x (p, q) 



p-q 




e 2 (p) 



Figure 11: The one- loop contribution to the self-energy the one-loop vertex correction 

Vi(p;q), and the two- loop contribution to the self-energy, ^{p). The calculation of £2(2?) 
can be broken down into two steps: first calculate Vi(p; q) and then replace one of the two 
bare vertices Vq which appears in Ei(p) by V\. 
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This odd part of the vertex correction has the form: 



Vi(p,e;q,^) = -c sgn(e) 



2 

e 2 +y 

OJ 



(8.59) 



where c is a constant; this behavior is found in the limit \lu\ > |e| Numerical calculations |94 
confirm the existence of this feature of the vertex, tracking the fermion energy |e| as pre- 
dicted. Linearization of the quasiparticle spectrum does not change the feature qualitatively. 
Naively equation p.59| would appear to be a subleading correction, but, because it has a dif- 
ferent symmetry than the bare vertex (which has no frequency dependence), KopietzpTf 
finds that it makes the leading contribution to the 2-loop self-energy, at least when bare 
fermion propagators Go are used, equation B.49 . If this is correct, there is no obvious sense 
in which a Migdal-type theorem holds. On the other hand, it would be surprising if Fermi 
liquid behavior were restored, as claimed by Kopietz and Castilla|98|]. In particular, it seems 
unlikely that the ln(e) term at two-loop order resums into a compensating power law as 



suggested in Reference |91| . We know that Fermi liquid behavior is not restored in the case 



of a linearized fermion spectrum, where multidimensional bosonization can be applied with 



confidence [92 



The form of the composite quasiparticle propagator therefore remains controversial. In 
any case, as the propagator is not gauge- invariant, it is unclear whether or not physical 
meaning can be ascribed to it. We turn now to an examination of some gauge- invariant 
properties. 



9 Two-Particle Response of Gauge- Fermion Systems 

In the previous section we focused attention exclusively on the single particle properties of 
the composite fermions of the Chern-Simons gauge theory. Using multidimensional bosoniza- 
tion we derived the composite fermion propagator which was found to have non- Fermi liquid 
form. However the composite fermion Green's function is not gauge invariant, and does not 



have a clear physical meaning. In fact it has been shown by Kim et a/.|)3| that, depending 
on the choice of gauge, different types of singularity appear. These same authors also used 
the perturbative loop expansion to study such gauge invariant quantities as the density and 
current response at small momentum transfer. These two particle response functions are 
gauge-invariant because the physical electron density, which is of course gauge-invariant, 
equals the composite fermion density. Kim et al. were able to show that the two-particle 
properties of the interacting gauge system are in fact Fermi liquid like: the singular correc- 
tions to the quasiparticle self-energy were canceled by singular vertex corrections, at least up 
to two-loop order. However the possibility of singular behavior at large momentum transfer 



could not be ruled out. In fact, Altshuler, Ioffe and Millis|9(J have argued that the vertex 
correction at momentum transfer 2kp should be singular since the gauge fields mediate an 
attractive interaction between currents moving in the same direction. In this section we use 
the multidimensional bosonization formalism to determine the density and current response 
for both small momentum transfer |q| w and |q| « 2kp. 



The density- density correlation function is given in all generality by equation |7.2| , in 
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particular at long wavelengths the density response is given by 

n 00 (q, u) = - £( J(S; q, u) J(T; -q, -w)) 

S,T 



(9.1) 



where in complete analogy to equation |7.11| the correlation function for the gauge theory can 
be written as 

(J(S;q,w)J(T;-q, -u)) 
n s • q 



+ 



A 



XD-1 

^ - v* F n s ■ q ' 

(n s • q)(n T • q) 



(2ir) 2 (iuj — Vph s ■ q}(iu — Vph T ■ q 



■b:(S-q)D^(q)b,(T-q) 



(9.2) 



Recall that D is the gauge field propagator, the vector b v (S; q) = (1, ii?^ q j^ s ), and the index 
v runs over (0, T). Now as 



6:(S;q)D^( g )6^(T;q) = A>o(g) - ^^^L>To(g) + iv* F ^^D 0T (q) 



iqi 

+ (q^ s )(qXn T ) 

|q| Iql 



q 



(9.3) 



only the term proportional to -Doo survives the summation over S and T and therefore 

II 00 (q, cu) = -xo(q) [1 - D 00 {q)xo(q)) (9.4) 

where 



A)o(<?) 



and 



detK = xo(q) 



Xr{q) - 
Xr(q) 



q^(q) 

(27T0) 2 . 

q 2 1/(q 



(detK)" 1 



q- 



(9.5) 



(9.6) 



(27T0) 2 J (2tT0) 2 

Here xo an d Xt are the longitudinal and transverse susceptibilities of the free electron gas 
defined in equations |6.48| and |8.20| . Combining these results, the long wavelength density 
response can be written as 



n o(q,^) 



n° o (q,u,)q 2 /(27T0) 5 



(9.7) 



q 2 /(27T0) 2 - ng (q, u)[n f /m* + g 2 V(q)/(27T0) 2 - U° n (q,u)} 

where we have introduced the paramagnetic component of the current response 

Xr(g) = n? 1 (g)-n / /m* (9.8) 

and I1q is the density response of the free electron gas. The transverse current response can 
be found in a similar way 



Il u (q, uj) = U° n (q, u) \l + D TT (q, u)U° u (q, u) 



(9.9) 
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where 



n§o(q,") 



qV(2vr0) 2 - U° 00 (q,u)[n f /m* + q 2 F(q)/(27T0) 2 - I^u;)] 



(9.10) 



The results found here for the longitudinal and transverse response, equations |9.7| and 



are identical to the results found within RPA||78||. They also agree with the leading two-loop 
result of Kim et al. from which we can conclude that the cancellation between divergent 
vertex corrections and divergent self-energy corrections is inherent in bosonization. The 
zero-frequency density response, the compressibility, vanishes as 



n oo (q,0) oc |q| 



(9.11) 



instead of |q| 2 or faster, so the electron system at v — 1/2 is in fact compressible, a re- 
sult reproduced recently within the formalism of Shankar and Murthy[[82| by Halperin and 
Stern (Tog. 

We now turn to consider the real space density response function in two dimensions at 
momentum transfer 2k p. With the gauge interaction, currents moving in the same direction 
are attracted to each other so we might anticipate nonanalytic behavior in n 2 fc F . This 
nonanalyticity has a quite different physical origin than the Kohn anomaly. The Kohn 
anomaly is due to the reduced availability of phase space at low energies. By contrast 
the behavior we expect in the present context should be apparent even if we only consider 
scattering between two patches at exactly opposite points on the Fermi surface. Therefore 
we need only consider terms with kg = — kx in the 2kp response function, equation |7.13| . 



n 2 fc F (x,r) 



- G^) ? < exp 



Att 



x exp 



n 



(S;O) + 0(-S;O)) 



(9.12) 



which can be transformed using the generalized Baker-Hausdorff formula, equation [2.26| , into 

2 



n 



2k F 



/ Q \ z 4-TT 

- (y-J exp -— ( [<0(S; x)0(-S; 0)) - (0(S; O)0(-S; 0))] 



- [(0(S;x)0(S;O)> 



(S;0))] + [S^-S] 



(9.13) 



The correlation function connecting patches at opposite points on the Fermi surface 
contains the new physics; it can be expressed in terms of the canonical boson correlation 
functions as: 



[(0(S;x)0(-S;O))-(0(S;O)0(-S;O))] 



n 

Air 



E 



duo 



1 



x 



P ,n s .p>cr 27r |ns-p| 
y(p.™-) _ i)( ( S ;p)o(-S; -p)> + (e^ p ™) - l)(a*(S;p)o*(-S; -p)> 



(9.14) 
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The boson correlation functions are then determined from the generating functional equation 

Q|n s • p| 



8.15 



(a*(S;p)a*(-S;-p)> 



(v* F h s ■ p) 2 - (iu) 2 



, 2 / p x n s \ 

Vp \ IpI J TT ^ P ' 



DM + iv F ^^-(D 0T (p) + £> T0 (p)) 



IPI 



(9.15) 



The other correlation function, (a(S;p)a(— S; — p)), is given by this same expression upon 
replacing iv* F — > —iv* F . Combining these results in equation |9.14| and adding the correspond- 
ing terms with S — > — S we find that the contribution to the exponent of equation |9.13| from 
the interpatch correlations, (S, — S), is given by 

~ {[(0(S; x)0(-S; 0)) - (0(S, O)0(-S, 0))] + [S - -S]} 



2E 



dw [e^ p - 



x— ut) 



1] 



2tt cj 2 + (i£n s • p) 



-D o(p,u) 



D TT (p,u) 



(9.16) 



The contribution of the intrapatch correlations can be obtained directly from the boson 
propagator, equation |3.21| : 



Att 

— MS; x)0(S; 0)) - (</> 2 (S; 0))] + [S - -S] = In ■ 



_ V f — [e i(p - x_LJr) - 1 



+ 7- 



n s ■ x + u^r)(n s ■ x - iv* F r) 
1 



x 



D„o(p^) + (^) : 



(iu - f^n s • p) 2 (iu + uf.ns • p) 2 



D TT (p,u) 



(9.17) 



and therefore 
n 2 fc F (x,r) 



A 



(27T) 



E'- 



1 



exp^E/^t^-^-l] 



x 



n s • x + if^r)(n s • x - iv^r) 
VFPjj-Doo(p) + 



' 1 IpI 



D TT (p) 



[(iu) 2 - (v F h s • p) 



212 



9.18) 



where the prime on the sum over patches indicates that the sum over S is only over patches for 
which |x x fig | A < 1. As in the case of the two-point function the most singular contribution 
comes from the p-limit, the exponent can be evaluated and we find 



n 2 fc F (x,r) 



A 



(27T) 



E' 



exp£(n s • x,r) 



J s 



n s ■ x + iv F T)(h s ■ x - ivpr) 



oc 



exp £( x 


,r) 


(|x 


+ iv F r)(\x 


— ivpr) 



(9.19) 
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where 



S(x,t) = C 
+ 



X 



X . ; v In | (x - m£t)<J| + t — - ; r 



In |(x + iv* F r)8\ 



^ In 2 | (a; - ivpT)S\ + - In 2 | (x + iv F r)5\ 



(9.20) 



for \x — w f t\5 ^> 1 and \x + iupr|5 1. The log-squared terms appearing in the exponent 
£(|x|, r) have their origin in the correlations induced by the gauge interaction between bosons 
in opposite hemispheres. Although it is technically difficult to Fourier transform the response 
function into frequency space, the leading singular behavior can be estimated by a rough 
scaling argument, and it is found to be 



U 2 k F {uj) oc exp[(ln 2 \u/v* F S\] 

where again v F S = O(ep). n 2 fc F is strongly singular as lj 
results from an exact diagonalization study JIT)]] ]. 



(9.21) 

0, in qualitative agreement with 



10 Nested Fermi Surfaces 

Understanding the cuprate superconductors, both in the normal state and in the supercon- 
ducting state, is one of the most challenging problems in condensed matter physics. Many 
interesting phenomena emerge, especially when the hole doping is less than optimal. As the 
doping is decreased superconductivity is quenched and eventually the material becomes an 
antiferromagnetic insulator; hence in the underdoped regime the interplay between the anti- 
ferromagnetic exchange channel and the superconducting channel is likely to be important. 
This is the common operating assumption of spin fluctuation [[102[], gauge || ||, |103|1 , SO(5) 



symmetry [|104|1 , and stripe |T05| models of high-temperature superconductivity, to mention 



just some of the theories. Indeed since angle-resolved photoemission experiments ||106|| indi- 
cate that the Fermi surface of the copper-oxygen plane is nested or nearly so; it seems clear 
that the tendency towards antiferromagnetism must be taken into account. 

The interplay between nesting, umklapp processes, and the formation of charge-density- 
wave (CDW), spin-density- wave (SDW), and BCS instabilities has been illustrated in a sim- 
ple model of the Fermi surface which reduces the nested planes to just four Fermi points^]. 
In a recent paper, Furukawa and Rice use a generalized version of this model to argue that 
the spin gap in the cuprate superconductors can be understood as a consequence of umk- 



lapp processes [TU7]. In this section, however, we retain the whole Fermi surface. Wilson's 
RG algorithm, which we introduced in Section |3], integrates out high-energy degrees of free- 
dom to derive a low-energy effective theory. In degenerate fermion systems this means that 
we shrink successively the inner and outer sides of the momentum shell which encloses the 
Fermi surface |108| |. As in Section ^ we investigate the RG flow of the effective interac- 
tions by integrating out the high-energy modes of the bosonic field <f>(S',x). The advantage 
of the bosonized RG calculation is that the existence of well-defined quasiparticles is not 
assumed a priori. Furthermore, interaction channels which possess {7(1)£° symmetry can be 



diagonalized exactly at the outset [30 
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First we introduce the prototypical low-energy model of a two-dimensional nested fermion 
liquid. Then multidimensional bosonization is used to incorporate the direct and exchange 
channels and to show that the resulting low-energy fixed point is a marginal Fermi liquid, 
which has obvious implications for the normal state of the cuprate superconductors [jS9fl . 
Next we derive the RG equations for the BCS channel and discuss the BCS instability. 
An analog of the Landau theorem, here made applicable to anisotropic Fermi surfaces, is 
invoked. We include the marginally relevant antiferromagnetic interactions and evaluate their 
RG flow and their contributions to the different BCS channels. The flow in the d^^-wave 
channel dominates the other BCS channels. This phenomena is similar to the Kohn-Luttinger 



effect fl47|, |109|| . The difference here being that nesting of the Fermi surface can induce BCS 
instabilities at leading order, rather than at subleading order as in the Kohn-Luttinger effect, 
and therefore can enhance the tendency toward superconductivity. 

The Hubbard model is a deceptively simple, but convenient, description of many strongly 
correlated fermion systems[ |110|| . It serves as the prototype for our analysis ||111| 1. A RG 
investigation of this model, of course, must take into account a realistic band structure near 
the Fermi surface. On a square lattice, the two-dimensional Hubbard model with nearest- 
neighbor hopping has a low-energy description which is complicated by several features 
of its anisotropic nested Fermi surface. Near half-filling, low-energy umklapp processes 
are permitted. The Fermi speed varies along the Fermi surface, and there are van Hove 
singularities in the density of states. Due to nesting, the system certainly cannot be a 
standard Landau Fermi liquid. The quasiparticles are expected to have a short lifetime, 
r ~ 1/T, and even the existence of a well-defined Fermi surface may be questioned |p.07| , |1 12 



Consider the Hamiltonian of the two-dimensional Hubbard model: 

H = -t 4%a + Uj2n x1 n xl , (10.1) 

<x,y> x 

where < x, y > are nearest neighbors on a square lattice and U is the on-site interaction 
strength. We assume \U/t\ <C 1 so that a weak-coupling perturbative expansion is justified. 
Of course in real systems typically \U/t\ ^> 1, but we assume as an operating hypothesis 
that the behavior at strong coupling is qualitatively the same as at weak coupling. This is 
equivalent to assuming that, at intermediate values of the coupling, no other fixed points 
intervene between the weak and strong coupling regimes. To begin, as usual, first we turn 
off the on-site interaction term and discuss the nature of the Fermi surface resulting from 
the hopping term alone. The spectrum of the fermions is 

e(k) = — It (cos k x a + cos k y a), (10-2) 

where a is the lattice spacing, not to be confused with 1/A, the ultraviolet cutoff of multidi- 
mensional bosonization. At half-filling, due to particle-hole symmetry, the Fermi surface is 
the set of points in k space with zero energy, and this set satisfies: 

cos k x a + cos k y a = . (10.3) 

The Fermi velocity depends on k and is given, on the Fermi surface, by 

vp(k) = 2ta (sin k x a, sin k y a) 

= v (sin k x a, sin k y a) (10-4) 
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where v = 2ta. Van Hove singularities in the density of states occur at the four corners of 
the Fermi surface, k = (±7r/a, 0) and k = (0, ±7r/a), as the Fermi velocity vanishes at those 
points. 

Now turn the interaction U back on and perform RG transformations until the low energy 
scale e c is reached. Assuming that the weak interaction does not alter the shape of the Fermi 
surface or change the spectrum qualitatively, we may linearize the quasiparticle spectrum: 

e(S;q) = v* \sm(k Sx a)\ q l{ {S) , (10.5) 

which agrees with equation pL0.2| close to the Fermi surface. Here q\\(S) = fis • (k — ks), where 



as usual ks is the Fermi wavevector for patch S and ks x is its x-component. This Fermi 
surface is difficult to bosonize since at the four corners, the van Hove points, the Fermi 
velocity is zero. However, away from half-filling the van Hove singularities are no longer 
in the low-energy zone. Therefore we introduce a simple spectrum in which the van Hove 
singularities are smoothed out. If we distort the four corners of the Fermi surface into four 
quarter-circles with a radius k as shown in figure [10), the modified spectrum is given by: 



v* | sin(k Sa .a)| <?||(S), k < \k Sx \ < vr - k, 
e(S;q) = <( (10.6) 
v* sin(Ka) ?||(S), otherwise 

This is a reasonable choice because the number of states near the van Hove points is of 
order na/ sin(fca) ~ 0(1), which is finite even upon taking the limit k — > 0. The smoothed 
quarter-circles contain a fixed fraction of the total number of states, as do the original van 
Hove points. From now on we set a = 1 and v* = 1 for convenience. Factors such as 1/a 
may be recovered whenever a comparison with other scales is needed. For example, with 
this convention the nesting vectors are Q = ±(7r,7r) and Q = ±(7r, — n) and the density of 
states at the Fermi energy, N*(0), which is of order (t>*a) _1 becomes a factor of order unity. 

Now that the Fermi velocity is well-defined over the entire Fermi surface, the model can 
be bosonized. As before, we coarse-grain the Fermi surface with small squat boxes of width 
A along the Fermi surface and height 

A(S) = -%r (10.7) 
v F {S) 

normal to the Fermi surface, as shown in figure [TD|. The Fermi speed in patch S is 

MS) = sin(k s .) . (10.8) 

The two momentum scales A and A(S) are taken to be small so that the inequality 1 ^> 
A 3> A(S) is obeyed in each patch. As we stressed in Section §, these limits are imposed 
to minimize the scattering of fermions outside of the squat box to conserve the number of 



fermions in each patch; the U(l) current algebra equation |2.11| holds in this limit. The 
effective action of the bosonized low-energy theory is 

S[a*, a] = J2 E / 7r[™ ~ MS)As • q] a* Q (S; q) a a (S; q) + Si[a*, a] . (10.9) 

q A(S)>n s -q>(r Z7T 
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Figure 12: The low-energy effective theory in the vicinity of the Fermi surface. The faces 
are flat, but the Fermi speed varies along each face. The corners have been rounded out into 
quarter-circles of radius k to smooth out the Van Hove singularities. The thin momentum 
shell is coarse-grained into squat boxes of height A(S) and width A. Again we require 
|q| < A(S) «A« l/a for bosonization to work. Momentum Q is one of four commensurate 
nesting vectors. 
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Here a a (S; q) and a* a (S; q) are the familiar canonical bosonic fields and Si is the interaction 
term. Now consider the renormalization of the fermion interactions. As in the case of a circu- 
lar Fermi surface, the marginal interactions are four-fermion (two-body) interactions, which 
can be verified by checking the engineering dimension. However, unlike the case of a circular 
Fermi surface, nesting permits additional marginal scattering channels at large momentum. 
In addition to the BCS instability, these channels could induce spin-density-wave (SDW) or 
charge- density- wave (CDW) instabilities; but first, before we discuss these logarithmically 
renormalized channels we must consider the forward scattering channel. We recall that, in 
the case of a circular Fermi surface, the interacting Fermi system with non-singular forward 
scattering is a Landau Fermi liquid in D > 1. In one dimension, however, forward scattering 
has a dramatic effect leading to a Luttinger liquid. Clearly the nested Fermi surface has 
some similarities to the one-dimensional system; it is important, therefore, to understand 
the effect of the forward scattering channel before the large- momentum scattering channels 
are included perturbatively. 

For simplicity, we suppress the momentum dependence of the coupling function and drop 
the spin indices. Inclusion of the spin exchange channel is not expected to change the result 
qualitatively. The action for the forward scattering channel in the fermion basis is 

S k \rM = -h E /^rfr^(S;x,r)V(S;x,r)^(T;x,r)^(T;x,r). (10.10) 

S,T J 



In bosonized form it is written compactly as 

fo ^ ^ r du 
s,t q 



* = W-*), 



(10.11) 



which is clearly marginal since it is bilinear in the boson fields. Here as usual q denotes 
the 2+1 dimensional vector (q, u). Focusing our attention on the fermions which lie on the 
nested part of the Fermi surface, since that is where the most singular behavior arises, we 
obtain the in-patch boson correlation function 



(a(S;q)a^(S;q)) 



1 



+ 



A 



n s • q 



iu — fi?(S)n s • q (27r) 2 [iu — i> F (S)ns • q) 2 
After analytic continuation to real frequencies, 

1 



[Mq)}' 1 • (10-12) 



K (q) 



2A ^ 0(ns • q) n s • q 



fo (2tt) 2 gU- v F (S)n s • q 

(2tt) 2 Xo(sin k q, u) + XN(q, u) 



1 V2 
+ 



fo (2tt) 2 

where Xn is the contribution from the flat part of the Fermi surface: 



(10.13) 



Xn(q) 



In 



1 — uj 2 jq\ + cosk 



^1 — uo 2 /q 2 i-yjl — u 2 /q 2 — cos K 
4 



In 



1 — uj 2 ' jq\ + cosk 



u 2 /q 2 



COS K 



10.14) 
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As before q» is the component of q parallel to the Fermi surface normal fig and q± is the 
perpendicular component. The contribution to the boson correlation function from the 



curved part of the Fermi surface, the term proportional to Xo i n equation |10.13| , as shown 
in Section || contributes a term to the fermion self-energy, Im E ~ uj 2 hiuj, and is therefore 
uninteresting. 

We can estimate perturbatively the imaginary part of the quasiparticle self-energy arising 
from the susceptibility Xn{q)- m the limits \q»\ > \u\ > \qn\ sinK or |g_J > |a>| > \q±\ smtc, 
the susceptibility has an imaginary part. The boson correlation function then can be ex- 
panded to second order in f : 

(a(S;g)a t (S;g)) 



u - v F (S)h s ■ q 



Just as for the isotropic Fermi liquid of Section [| the first-order term in fo renormalizes the 
Fermi velocity, and the velocity correction, 5v = /oA, is infinitesimal in the A — > limit. 
Spin-charge velocity separation therefore does not occur, as can be verified easily by inclusion 
of the spin indices. Inserting equation |10.15| into equation exponentiating, expanding 



the exponential to second order in /o, and finally performing a Fourier transform into (k, u) 
space gives the imaginary part of the quasiparticle self-energy|30| , Im £(k, uj) ~ uj /q A, which 



has the marginal Fermi liquid form] 89 1. Like the D = 3 case of a spherical Fermi surface, 
the prefactor A which appears in the self-energy obtained by bosonization is replaced by 
a number of order unity in an RPA calculation of the self-energy. Corrections at higher 
orders in f can be estimated and are also consistent with the marginal Fermi liquid form, in 
agreement with expectations of others who found that the inverse quasiparticle lifetime, 1/r, 



is proportional to the temperature [|112|| . Since the quasiparticle weight in a marginal Fermi 
liquid vanishes logarithmically at the Fermi surface, the Fermi liquid fixed point breaks down 
but the occupancy still changes rapidly enough near the Fermi surface for multidimensional 
bosonization to be applied. 

Now we include the large-momentum scattering channels by the perturbative RG method. 
We can discard the forward scattering channel since we have shown that it contributes only 
sub-leading corrections to the RG equations. We study the RG flow of the BCS channel 
and derive an analog of the Landau theorem [|113|| . Like the circular Fermi surface, the 
nested Fermi surface has reflection symmetry about the origin and therefore the BCS channel 
flows logarithmically with decreasing energy scale. For a circular Fermi surface, the Landau 
theorem states that there is a BCS instability if any one of the Ve is negative, where Vi is the 
£-th orbital angular momentum component of the BCS vertex function defined in equation 



3.13| . However, for nested Fermi surfaces, the BCS vertex is not rotationally invariant, and 



Vbcs(S, T) must be expanded in a double Fourier series |109| and the Fourier coefficients are 
now coupled. The RG equation for the BCS flows, equation [3.11| , is also modified as the 
Fermi speed Uf(U) is not uniform. 

dV BCS (S,T) A Vbcs{S,V) Vbcs(U,T) 

d\n(s) (2vr) 2 ^ v F (U) ' 1 ' ' 
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Note that the weight function ^(U) -1 now appears in the sum. 

In order to analyze the RG equation, we must Fourier-decompose the BCS vertex func- 
tion. The vertex function has -D4 symmetry, and there are five sets of basis functions which 
generate the two-dimensional function space to which it belongs. The five sets are: A\ and 
A 2 (s-wave), B\ and B 2 (d-wave), and E (p-wave) ||109|| . The basis functions for each set are 
as follows: 



A l : 


ai(n, S) = 


A 2 : 


a 2 (n, S) e 


B x : 


h(n, S) = 


B 2 : 


b 2 {n, S) = 


E : 


e(n, S) ee 



MS) 1 1/2 



2?r\/2 



MS) 



cos(2nk Sx ) , ai(0, S) 



v F {S) 



1/2 



MS) 



1/2 



sin(27ik Sa .) 



1/2 



[Msr 


1/2 ^ 


.2%y/2. 





cos[(2n + l)k s J , 
MS)! 1 ^ 

£v^\ sin[{2n + 1)ksi 

ci sin[(n + l/2)k Sx ] 



+ c 2 sgn(k Sy ) cos[(n + l/2)k Sa 



10.17) 



A factor of yvp(S) has been inserted to render the functions orthonormal under summation 
over the Fermi surface with the weight function uf(S) -1 . In general the BCS vertex function 
can be expressed as a double Fourier expansion in terms of these orthonormal basis functions; 
there are no cross terms between basis functions of different symmetry: 



BCS 



'S,T) 



]T {C? m oi(n, S) ai {m, T) + a 2 {n, S) a 2 {m, T) 



n,m 



+C 3 " m h(n, S) him, T) + b 2 (n, S) 6 2 (m, T) 
+C 5 " m e(n,S) e(m,T)} . 



(10.18) 



This expansion is valid as long as vp(S) 7^ which is the case for the model considered here 
as the Van Hove singularities have been smoothed out. 



Inserting equation |10.18| into equation |10.16| the RG equations for the Fourier coefficients 
are found: 



dC? m 
cHn(s) 



(2 



7T 



1 in 



(10.19) 



1=0 



where i = 1, 5 labels the 5 sets of coefficients. The analog of the Landau theorem for the 
nested fermion liquid is obtained by considering the diagonal Fourier coefficients which obey 
the equations 

The diagonal coefficients decrease as the length scale s increases, and therefore, if any one 
value of Cf n is negative initially, a BCS instability occurs. This sufficient condition for an 
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Figure 13: A typical density- wave scattering channel which is permitted because Qj + Qj is 
a reciprocal lattice wavevector. 

instability allows us to identify the possible symmetry of the resulting order parameter, since 
the gap function A(S) is related to the BCS vertex function as usual by: 

A (S) =-£ Vbcs{8, T) , (10.21) 

where E(\t) is the dispersion of the low-lying quasiparticle excitations in the superconducting 
state. 

With these criteria we can study the possibility of a BCS instability in the presence of a 
repulsive Hubbard-£7 interaction. First, we consider the RG flow of the channels generated 
directly by the repulsive-?/ interaction. We consider only those channels which flow logarith- 
mically in length scale s, as all other channels are irrelevant. In addition to the BCS channel, 
there are three other classes of interaction channel which renormalize logarithmically in the 
presence of nesting. The first class is the CDW or SDW channel, which we denote simply by 
DW (density wave). In the physically relevant case of a repulsive bare Hubbard interaction, 
the DW channel leads to a SDW instability. However, this channel is decoupled from the 
BCS channel. The second and third classes are the Q = and the umklapp Q ^ repulsive 
exchange channels which do renormalize the BCS channel, and which therefore can induce 
superconducting instabilities. 

Expressed in terms of the fermion fields, the action for the DW channel is: 

S,Ti=lJ=l 

lT(T; x) V*(T + Q,; x) [SPjP - 5%] . (10.22) 



x 



Here Qj and Qj are two of the four nesting vectors shown in figure [L^, i,j = 1,2,3,4, and 
the notation S + Qj is shorthand for momentum k s + Qj. Umklapp scattering, Qj + Qj ^ 
0, is permitted because Qj + Qj is a reciprocal lattice wavevector. The spin indices in 
equation |10.22| have been symmetrized consistent with Fermi statistics [IDS|. Due to nesting, 



this channel flows logarithmically with length scale and leads to either a CDW or a SDW 
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instability, depending on the sign of fnw- The action can be decoupled into CDW and SDW 
channels, 

If i 4 

Sdw = o / drd2x {fsDw S(S, Qi, x)-S(T, Qj, x)+f CDW n(S, Q;; x) n(T, Qf, x)} , 

S,T i=lj'=l 

(10.23) 

where 

S(S, Q i5 x) = i V* a (S; x) ^(S + C^; x) (10.24) 

and 

n(S, Q,; x) = V* a (S; x) ^ a (S + Q*; x) (10.25) 

are the spin density and the charge density, respectively. The primed summation over S, T 
indicates that the sum is only over the flat parts of the Fermi surface. 

Consider the RG flow of the CDW term; for simplicity again we drop the spin indices. 



From the bosonization formula, equation |2.21| , the bosonized form of the CDW action is 

c rn fcDW A 4 w 4 r 2 A(S) A(T) 

ScDw[4>\ = }^ L J drd x 



S,T i=lj=l' 



r \l All i 

x exp O^-[0(S + x) + 0(T + Q^; x) - 0(S; x) - 0(T; x)]} (10.26) 

Under rescaling, t — > t'/s, x\\ — > x',</s, and e c — >■ se' c ; as expected Scdw is scale invariant 
and therefore marginal. Next the fields are split into high and low-energy modes 0(S; x) — > 
0'(S; x) + /i(S; x) where, as in Section ^, /i(S; x) is the high-energy field with momentum in 
the range A'(S)/2 < |p| < A(S)/2 with A'(S) = €' c /vf(S). The correlation function of the 
h(S; x) fields is given by 

(/i(S;x,r) /^(T;0,0)) = 



S d-i £ ln ( , | fis x x|A « 1 



0, |n s xx|A>l, 

(10.27) 

similar to equation [3l| but with explicit patch-dependence in the speed ujr(S) and in the 
cutoff A(S). We then determine how the partition function renormalizes as the fast modes 
are integrated out: 

(exp [ - Scdw]) = (1 - Scdw + ^ S cdw + • • •) 

= e W [-S' CDW ) . (10.28) 

Here the average in the brackets is taken over only the high-energy fields h(S; x) and S' CDW 
is the renormalized CDW action. We evaluate (Scdw) using equation |3]^ and the formula 

(e ih ) = e~^ 2 / 2 > to obtain 

/q \ _ fcDW , A n2W >p f , ■* A(S) A(T) 

{ cDw) " ^ { Wy ] t Tl ^J & 

x exp {i^WiS + Q,; a) + 0'(T + Q,; z) - 0'(S; x) - <//(T; x)]}(10.29) 
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Replacing A(S) with sA'(S) [A'(S) is the new momentum cutoff] and 0' — » we recover 
equation |10.26| with the new cutoff. The third term on the right hand side of equation |10.28 
gives the second order correction to Sqdw- Using the relation (e A e B ) = e( AB+ ( A +B and 
equation |10.27| , we find 

A l 4 ^^f, , A(S) A(T) A 2 (U) 



fr 



CDW 



(27T) 



^2 J ' dr d 2 x 



S,T,U 



16s 4 



x exp < i 



{< 



n 



x 



i/(T;x)]} 



x 



J da d 2 y exp {^[(h(U; y, a) h(U; 0, 0)) + (h(U + Q,; y, a) h(U + Q ; ; 0, 



+ subleading terms 



2 

CDW 



A 



(2k) 



IV' ^ t J, J2 W A2 ( U ) 27r [ A A(~ 

J2 E drdx ^4^T dad{niJ ' y 

J S,T i=lj'=l U * 



4s 2 /A 2 (U) 



S,T i=lj=l 
47T 



K(U)a] 2 + (nu ■ y) 2 



A(S) A(T) ; 
+ subleading terms . 



i'(S + Qi, x) + 0'(T + Q i; x) - 0'(S; x) - 0'(T; x)\ } 



(10.30) 



Performing the integrals over fiu • y and a, we obtain the RG equation for the CDW coupling 
function 



dfcDw(s) 



A 



1 



fl 



(10.31) 



dln(s) 2(2tt) 2 ^ v f (U) JCDwy 

This channel is unstable if fcDW is negative, which is the case for a bare attractive interaction, 
U < 0, leading to charge-density ordering at wavevector (k,k) or (k, —k). 

Similarly, a RG transformation can be performed to find the flow equation in the SDW 
channel yielding 



dfsDw(s) 



A 



1 



fsDw( S ) ■ 



10.32) 



tZln(s) (2tt) 2 ^ ^(U) 

There is an instability if fsDW < 0, which corresponds to a bare repulsive interaction U > 0. 
The SDW instability leads to antiferromagnetic spin-density ordering, also with wavevectors 
(k, k) and (k, —it) as the repulsive interaction induces an antiferromagnetic exchange interac- 
tion. We can determine easily the energy scale of the SDW instability. Integrating equation 



10.32| , we obtain the critical energy at which the SDW instability occurs by equating it to 

(10.33) 



the energy scale at which the coupling constant fsDW grows to be of order one: 

1 



E S dw ~ E F exp 



(2tt)2 



ln f I/' 



SDW I 



Now we can examine the overlap between the BCS and the DW channels. In figure [14 
we show a typical scattering process belonging to both channels. The contribution of these 
processes to the action can be written in terms of the BCS vertex function and the fermion 
fields as: 



Overlap = E/ rfr ^ VbCs(S, T) ^(Sj x) ^(S;x) M~ T ^) M^x) 4.S+Q ■ 

(10.34) 



S,T ' 
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Figure 14: A typical overlap between density-wave and BCS channels. 




Figure 15: (a) An exchange channel with zero momentum transfer which logarithmically 
renormalizes the BCS channel, (b) Same, but with momentum transfer 2Q. This is an 
umklapp process. 

Clearly, due to the factor of 5t,s+q; there is a severe phase space restriction in the overlapping 
channel, which gives rise to an overall factor of Aa <C 1. Therefore, at leading order the DW 
channel is decoupled from the BCS channel. 

The two classes of repulsive channels which do renormalize the BCS channel are shown 
in figure Their contributions to the action, written in terms of the fermion variables, are 
respectively: 

S R [r a ,M = |EEE/^ A r(S!x)^(S + Ap+Q !; ,) 

Z S,T i=l Ap J 
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x^* 7 (T;x) ^(T-Ap-Q^z) 
which is a zero momentum process and 



8 

a"7 



(10.35) 



E E E / dr d 2 x r a (S; x) M$ + Ap + Q, ; x) 



fu 
2 

^ S,T i=l Ap ' 

x^ 7 (T;a;) ^(T - Ap + Q-; x) 



a "7 ^y^a] ) 



10.36) 



which is an umklapp process permitted as 2Qj is a commensurate wavevector. In both chan- 
nels we assume the coupling functions to be momentum-independent and we anti-symmetrize 
the spin indices. Now we expand the partition function to second order in Sr + Sjj, inte- 
grate out the high-energy modes, and study the resulting effective action generated by this 
procedure. In the course of this RG procedure, both Sr and Sjj are renormalized as well as 
the BCS action, and the resulting flow equations are a complicated set of coupled differen- 
tial equations. As we are only interested in the initial growth of the BCS channel, we do 
not need to perform a complete analysis of the flow equations, we need only examine how 
Sr + Sjj renormalizes the BCS channel in the initial stages of its growth near s = 1. The 
second-order correction to the BCS channel is given by: 



1 



((Sr + Sjj) 2 ) 



Ms) (/£ + 4#)E£E/dr d 2 x r a (S;x) ^(-S;x) 

S i=l Ap J 

x^(S + Ap + Q,; x) ?/U-S - Ap - Q i; x) 
\/2 1 , ( \ + cos(k + |Ap x |/2)- 



2(2tt) 2 cosAp x ./2 



In 



, ln ( g ) j2 



R 



V2 



(27T) 



hi 



1 — cos(k 

1 + COS K * 
1 — COS K 



|Apx|/2) 



1 

2^ 



x 



Y,j dr d 2 x V* a (S; x) ip*^{S; x) ^(-T; x) ^ Q (T; x) 



S,T 

+ non-BCS terms . 



(10.37) 



The second term in equation |1 0. 3 7| governs the flow of the s-wave BCS channel but since the 
s-wave channel is repulsive at first order it is not a candidate for instability. However, the 
first term in equation |10.37| , as it has strong wavevector dependence, contributes to the RG 
flow in all of the other BCS channels. The flow equations are: 



dV, 



BCS 



'S,T) 



dln(s) 



A 



V2 



UJ 2(2ir) 2 cosA P;E (S,T)/2 



In 



1 + cos(k + |Ap x (S,T)|/2) 



E 



(2vr)2 ^ Vf (V) 



l-cos(«+|Ap :B (S,T)|/2)' 
V BCS (S,U) V BCS (U,T) , (10.38) 



if S and T are in the same quadrant in the (k x , k y ) plane. On the other hand, when S and 
T are in different quadrants, 



dV, 



BCS 



:s,t) 



A 



dln(s) 



E 



V B cs(S,U) V BCS (U,T) 



[10.39) 
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Finally, we can determine which BCS channels are attractive by evaluating the generalized 
Fourier coefficients of the double expansion shown in equation |10.18| . After performing 



the generalized Fourier transform of the RHS of equation |10.38| we find that, among the 
fundamental modes, there are two diagonal Fourier coefficients which flow to negative values. 
Their rate of change, evaluated for the case k = 0.01, is: 



dln(s) 



dCf 
dln(s) 



1 



^ MS) v F (T) 



dV, 



BCS 



S,T) 



0-1 (/I + 4/, 
1 



a 2 E 



dln(s) 

^T/(2tt) 2 . 
dV BCS (S,T) 



^ v F (S) v F (T) dln(s) 
0.23 (/ 2 +4/ 2 )-|C 3 00 | 2 /(27r) 2 



a 2 (l,S) a 2 (l,T) 



6i(0,S) 6x(0,T) 



(10.40) 



(10.41) 



Here we have set C\ 1 {s = 1) = C®°(s = 1) = initially. From the Landau theorem we infer 
that there can be a BCS instability in either the B\ (d x 2_ y 2-wave) or A<i (extended s-wave) 
channels provided the SDW instability does not preempt it. Furthermore, as the inequality 



dCf 



dln(s) 



> 



dCl 1 



dln(s) 



(10.42) 



holds for small values of k, d x 2_ y 2-wave superconductivity is preferred. Superconducting 
order, if it exists, will most likely possess 0^2 _ y 2 -wave symmetry. This conclusion is consistent 
with the expectation that nearest-neighbor antiferromagnetic exchange favors d x 2_ y 2-wave 



superconductivity. It also agrees qualitatively with a diagrammatic approach [ |114| | and a 



parquet solution obtained for the flat Fermi surface [115 



Of course, the tendency towards spin-density order dominates all other instabilities close 
to half-filling. But slightly below half-filling, because nesting is not perfect but nonetheless 
effective, there is competition between the SDW and BCS instabilities. As the energy cutoff 
e c around the Fermi surface is reduced via the RG transformation towards e s , where e s is 
the energy deviation of the actual Fermi surface from perfect nesting, the SDW channel 
stops flowing but the BCS channel continues to renormalize logarithmically as its flow is 
independent of nesting. If the initial SDW coupling is small enough or if e s is large enough 
(the SDW channel does not develop an instability until e c — > e s ) SDW order will not occur. 
Instead, superconductivity will set in at a sufficiently low temperature. 



11 Conclusion 

In this review we have shown that multidimensional bosonization, when combined with 
the Renormalization Group picture, is a powerful tool for the study of strongly interacting 
fermion systems. First, an effective action for quasiparticles is derived by integrating out the 
high energy degrees of freedom which lie outside of the concentric shells which encircle the 
Fermi surface. The renormalization process is stopped before all of the degrees of freedom 
are eliminated; the effective Hamiltonian so obtained describes quasiparticles confined to a 
narrow region around the Fermi surface of width A <C kp. Of course no attempt is made to 



78 



calculate the values of parameters such as the effective mass m* or the coupling constants 
for the short-ranged part of the forward and exchange interactions as this can only be done 
reliably in special limits, for example the case of a very weak bare interaction. The Fermi 
surface is now discretized, tiled into squat boxes with aspect ratio A < A, where A is the 
length of the box along the surface. In each squat box, currents defined in terms of the 
fermionic quasiparticles also have an equivalent representation in terms of bosonic fields. 
When it is legitimate to linearize the quasiparticle spectrum about the Fermi momentum by 
dropping the non-linear terms, the complicated many-body interactions reduce to a quadratic 
form in the boson fields which possesses the large U(l)f symmetry, assuming that there is 
no BCS instability. Given these caveats, the boson propagator can be determined exactly 
from a generating functional; the fermion propagator then follows on using the exponential 
representation of the fermions in terms of the bosons. 

As the existence of Landau quasiparticles is not presupposed, multidimensional bosoniza- 
tion is a useful nonperturbative approach to access non-trivial ground states of strongly cor- 
related fermion systems and to study their stability. In one spatial dimension, well-known 
exact results are recovered. The ground state is generically a Luttinger liquid, and the 
low-lying excitations are collective modes of charge and spin which propagate at different 
velocities. In higher dimension the Landau Fermi liquid is realized when fermions interact 
via a longitudinal interaction of either short-range, or of the usual Coulomb form. Unphysi- 
cally longer ranged interactions, however, can destroy the Landau quasiparticles. All of the 
important results of Fermi liquid theory can be recovered within bosonization, the nonana- 
lytic T 3 In T contribution to the specific heat in three dimensions, the equations of motion 
of charge and spin collective modes, as well as screening and quasiparticle damping in two 
and three dimensions. Both long and short wavelength physics can be described provided 
the energy is low enough to be within a shell surrounding the Fermi surface. For example 
the density response has the usual RPA form at low wavevectors, and at high wavevectors 
near 2k p the Kohn nonanalyticity is found. 

The possibility of non-Fermi liquid ground states has been raised for systems interact- 
ing via transverse gauge fields, as these fields are not screened at low frequencies. Such a 
scenario has been suggested as the mechanism underlying cuprate superconductivity, but 
a concrete example is furnished by quantum Hall systems in two dimensions, specifically 
a two-dimensional electron gas in a perpendicular magnetic field such that only the lowest 
Landau level is half-filled. In the widely used Chern-Simons gauge theory the physical elec- 
trons are mapped by a singular gauge transformation onto composite fermions which each 
carry two units of statistical flux and interact via a transverse gauge field in addition to 
the usual longitudinal density-density interaction. At the mean field level the average flux 
cancels the external magnetic field and the ground state is a free Fermi gas. The stability 
of this state against gauge (= charge) fluctuations was studied by determining the com- 
posite quasiparticle propagator. Multidimensional bosonization broke down for longitudinal 
interactions shorter-ranged than Coulomb as the transverse gauge fluctuations, which are 
moderated by the longitudinal interaction, are too strong in this instance. In the case of the 
Coulomb interaction, two important results were found. At equal time and at long distances 
the propagator has the form Gjr(S; xn) ~ In the opposite limit of long temporal 

separation, \t\ ^> the leading terms can be Fourier transformed and the self-energy 
extracted, yielding the marginal Fermi liquid form, T,(u) ~ C^ 11 M> identical to that found 
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within RPA. 

This result has been criticized by several authors and remains controversial. The basis 
for the criticism is the linearization of the quasiparticle dispersion, an approximation which 
is especially useful in bosonization as it leads to a Gaussian bosonic action. Because gauge 
fluctuations carry relatively large transverse momentum, it has been argued that non-linear 
corrections to the quasiparticle dispersion should be retained, especially in the case of short- 
range longitudinal interactions. Furthermore, Altshuler, Ioffe, and Millis|5lJ argue that, in 
the case of the Coulomb interaction, when the quadratic q\j (2m*) term is retained the MFL 
form of the self-energy given above is the leading term, and all other higher-order terms 
are subleading. This claim has been criticized by Kopietz who points out|9T| that the two- 
loop self-energy is more divergent than the single-loop self-energy. Kopietz and Castillap8 



attempted to incorporate the q]_/2m* term within a bosonized treatment and find that Fermi 
liquid behavior is recovered, even in the case of short-range longitudinal interactions where 
the gauge fluctuations should be strongest and most disruptive. The nature of the composite 
quasiparticle propagator, and whether or not it is even physically meaningful, remains the 
subject of ongoing research. Different, direct, schemes for bosonization of Landau levels have 
been proposed ||116| , 117 |, and these may prove useful. Renormalization group methods, when 
combined with the Ward identity approach [ ITS | are another route being pursued. 

A particularly important feature of multidimensional bosonization is that it can be com- 
bined with the renormalization group to analyze various instabilities of the Fermi liquid, 
all of which are treated on an equal footing. Bosonization is most useful in this context 
as the forward-scattering parts of the interaction, those which preserve the U(l)^ symme- 
try, can be incorporated exactly at the outset as part of the Gaussian fixed point. As a 
simple example we have studied the BCS instability for a circular Fermi surface, and have 
rederived the RG flow equations found by Shankar in the fermion basis. For a nested 
Fermi surface close to half-filling, we studied the competition (for a repulsive bare interac- 
tion) between density wave and superconducting instabilities. A pleasing conclusion coming 
from this calculation is that the d x 2_ y 2 superconducting channel is favored in agreement with 
experiments [TT9|]. Above T c the normal state is a marginal Fermi liquid, also in qualitative 
agreement with many experiments ||120[| . Bosonization is particularly appropriate here as the 



forward-scattering channel responsible for the marginal Fermi liquid state is incorporated 
non-perturbatively. Unlike other approaches, in multidimensional bosonization Feynman 
diagrams receive rather less emphasis than algebraic and geometric ways of understanding 
correlated fermions. When combined with the renormalization group, bosonization yields in- 
sight into fermion liquids which does not rely on the resummation techniques of microscopic 
Fermi liquid theory P3|. By unifying the microscopic and Landau pictures, multidimensional 
bosonization provides a new way to think about correlated fermions. 

Finally we should mention some outstanding problems which remain unanswered. First, 
as we have noted, linearization of the fermion dispersion is a central approximation, crucial 



to the success of the approach. As pointed out originally by Haldane||i7||, nonlinear cor- 



rections to the dispersion can only be included systematically in a perturbative treatment. 
When these terms are important, it is not clear how useful bosonization remains. The most 
important failing of bosonization to date, however, is that it remains unable to account for 
the intermingling effects of interaction and disorder leading to diffusion and the modification 
of screening ||1 2 To understand the role of disorder it is necessary to confront the effects 
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of interpatch scattering, which breaks the U(l)f? symmetry down to just the usual global 
£7(1) charge conservation symmetry. Some progress has been made in the special case of 



one spatial dimension where the weak-coupling RG equations have been worked out||122 



In higher dimension, forward scattering by impurities can be included rather easily |38|, |123 



but this has no effect once the X/kp — > limit is taken, as required if multidimensional 
bosonization is to be rigorous. In the regime of strong disorder, Belitz and Kirkpatrick 
have used a Hubbard-Stratonovich transformation to decouple the fermion interaction along 
with replicated disorder by introducing bosonic Q- fields [ [124|| . They emphasize that the soft, 
low-energy, modes of disordered systems differ qualitatively from the low-energy degrees of 
freedom in the clean limit, as they are diffusive rather than ballistic in character. It remains 
to be seen whether or not a theory can be constructed which interpolates accurately between 
the limits of weak and strong disorder. 
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A Calculation of the Boson Correlation Function 



Use equation 2.18 to find the equal-time correlation function 



GV(S;x) = (0(S;x)0(S;O)-0 2 (S;O)) 



Air 



E ^(S;p: 



p,n s p>0 



/gip-x _ ]^ e -an s -p 

|n s ■ p| 



(a.i; 



where we have introduced the convergence factor e _ans p . To evaluate the integral over the 
component of momentum parallel to the surface normal, p» = ns • p, consider first the 
integral 

rV2 dp\\ e iv w x \\ 



here a is an infrared cutoff which drops out of the calculation at the end. Then 



'<o) = s i» 



2^ 



dl(xu 
dxu 



X/2 ^Pl e ip n (x }l +ia) 

4ir 



4ir(x\t + ia) 



(A.2) 

(A.3) 
(A.4) 
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and 



I(x\\ 



where, from equation |A.3|, 



— lnfxii + id) + C 

4tt 11 



4% 2a 



and therefore 



Kxn 



Air 



In 



2a x\\ + ia 
A ia 



(A.5) 

(A.6) 
(A.7) 



Now consider the integrals over the momenta perpendicular to the surface normal, special- 
izing to the case D = 2 for concreteness: 



A / 2 dp A 



ip±x A 



-A/2 2tt 



1 . Ax ± 
sin 



71X± 



A|x ± |«l 



, A\x±\ > 1 



(A.i 



Combining these results, and using the fact that J2 P — L D J we find that the infrared 
cutoff a drops out and 



GJS: x) 



"fin 



— oo, 



xu+ia 



Alx, I < 1 



A|x ± | > 1. 



(A.9) 



B Integrating Out Fast Boson Fields 

To evaluate the propagator of the h fields, make use of the relationship between <fi and the 



canonical bose fields, a and a*, equation 2.18 to write 



(h(s,x)h(s,o)) = £ f^'K^)^^))^) 



q\\ 



where 



and 



^ Gtt) I-A/2 d U L 
(o(S,g)o*(S,?)> 



A/2s<|q N |<A/2 
-1 



2cj - v F q\\ 



(B.l) 

(B.2) 
(B.3) 



The frequency integral is straightforward and equation |B.1| reduces to 

(h(s,x)h(s,o)) = — d n [o(n)e(r) - e(- n )e(-T)) (ba) 

47T J \/2s<\ qil \<\/2 q\\ 
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where we have used the result (given here for D = 2): 

*% ±e ^ = fl^J A ' AM<<1 (B.5) 
A ^ 2 l x -L/ ) q ; otherwise 

To evaluate the integral over gn write it as the difference J A ^ 2 — J A//2s an d replace the upper 
limit by a smooth cut-off to find 



(h(S,x)h(S,0)) 



fin 



X||+raJ,r+(2is/A)sgn(T) 
a; || +raJ,r+(2i/A)sgn(r) 



A|x ± | < 1 

(B.6) 

otherwise . 



C Free Fermion Propagator 



Expressed in terms of the boson fields the fermion propagator in patch S is given by 



G F (S;x, r) 



-A 



D-l ^tks-x 



(2tt) d a 



exp 



Air 

— <0(S;x,r)0(S,O)-0 2 (S,O)> 



(c.i; 



where the correlation function of the <fi fields expressed in terms of the canonical boson fields 
a, and a*, is given by 



G>(S;x, r) 



(0(S;x,r)0(S,O ) O)-0 2 (S;O,O)) 



47T 



E 



2^ 
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n s ■ q 



(C.2) 



which determines the fermion Green's function once the Fourier transform of the boson 



Green's function is known. Now the free boson Green's function is given by equation p.3| : 

-1 



and therefore 



G^(S;q,w) 



in 2 /-A/a 



iuj - Vpiis ■ q 



Gj(S;x,r) = — / dg. 



(C.3) 
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where we have used equation B.5 , valid for D = 2, to integrate over the perpendicular 
momentum. Now consider 
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and therefore 
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but lim r ^ + G^S; 0, r) = 0, which gives C = ^ ln(ia) and 



G°(S;x,r) 



47T 



In 



XII + ZUpT 



/(-/ 



For the case of D = 2, it follows from equation |C.1| that 



G F (S;x,r) 



iA 



(27r) 2 n s • x + WpT 
The Green's function for real time is obtained by continuing it 



(C.7) 



D How Kohn's Theorem Is Satisfied 

In this appendix, we derive the result quoted in Section g, namely the collective excitations 
of the interacting quasiparticles of the half-filled Landau level occur at a frequency uj = 
27rrif<j)/m, the cyclotron frequency of a particle with bare mass m in a magnetic field B = 
27rn/0/e. To prove this result, we add to the action of the Chern-Simons gauge theory a 
Fermi liquid interaction in which all the coefficients except f\ have been set equal to zero [[81] . 



The contribution to the action, equation |8.13| , is given in real time by 

fx 



2,kp 



d 2 x dt ip*(x)Vip(x) ■ ^*(x)V^(x) 



(D.l) 



which is made gauge covariant by replacing V by V — iA. On Fourier transforming it is 
given in terms of the currents as 



S 



h 
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2Vk 
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F q 



As in Section || the term quadratic in J is decoupled by introducing additional transverse 
A t and longitudinal Ai gauge fields. 

Including the Fermi liquid interaction the complete action of the theory is given by 

S[A",a,Z,e\ = E E /^{(«-i;Jfls.q)o*(S; ff )o(S;g) 

S q,nsq>0 

+ e(S;?)o*(S;g)+r(S;g)o(S;g)} 

M-i) + ^ 4Hi>} 



S q 

n s ■ q 
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- l^jAh)-}^)^-?)}, (D.3) 

where /i runs over 0,T,l and t. We have coupled the boson fields a* and a to £ and £* to 
construct a generating functional as in equation ^.15| . After integrating out the boson fields 
we obtain the effective action: 

+ 2 M A (q) A T {-q) + [-1 + xi(q)} Mq) A(-q) + [-i + Xt (g)] A t ( 9 ) A t (-g) 

27T0 jl jl 

+ xro(g) [4(g) 4(-g) - 4(g) + 2[v* F + f -f^} xM) A T {q) A t (-q) 



+ Xo( q )Mq) M -q)}-± £ /^ E{ ^ S ^V^s-ql 
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and the susceptibilities Xa(q) are defined by: 
/ s , rtM / cos0 
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and x = u>/(vp\q\). The boson propagator is obtained in the usual way 

(a(S;q) a f (S;g)) = — — 

u; — t>^n s • q + ir] sgn(u;) 



x 
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and the collective modes are given by the poles of the gauge propagator. The inverse gauge 
propagator can be read off from the action equation [D.4j and in the limit of interest uo ^> v* F \q\ 
is given by 



[K{q)\ 



/ N*(0) 
' 2x 2 

ilsL 

2x 
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\ 



) 



:d.s) 



The poles of D a/3 are found by setting det|ii" a( g| = which gives u 2 = [27r0riy(l/m* + 
/i/47r)] 2 = (27rrif(j)/mb) 2 - The second equality follows from the application of Galilean 
invariance to a Fermi liquid|f29[. The energy of the collective mode is at the cyclotron 



frequency determined by the bare electron band mass in agreement with Kohn's theorem. 
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